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Abstract 

We review the present state and future outlook of our understanding of neutrino masses and 
mixings. We discuss what we think are the most important perspectives on the plausible and 
natural scenarios for neutrinos and attempt to throw light onto the flavor problem of quarks and 
leptons. This review focuses on the seesaw mechanism, which fits into a big picture of particle 
physics such as supersymmetry and grand unification providing a unified approach to the flavor 
problem of quarks and leptons. We argue that in combination with family symmetries, this may 
be at the heart of a unified understanding of the flavor puzzle. We also discuss other new physics 
ideas such as neutrinos in models with extra dimensions and possible theoretical implications of 
sterile neutrinos. We outline some tests for the various schemes. 



* With permission from the Annual Review of Nuclear and Particle Science. Final version of this material 
is scheduled to appear in the Annual Review of Nuclear and Particle Science Vol. 56, to be published in 
November 2006 by Annual Reviews (www.annualreviews.org). 
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I. INTRODUCTION 

The Standard Model (SM) of electroweak (EW) and strong interactions is an overwhelm- 
ingly successful theory for particles and forces. It probes physics below a 100 GeV's (in some 
cases up to several TeVs ) and has met the challenge of many high-precision experiments. 
There are, however, strong reasons from considerations of both particle physics as well as 
cosmology to suspect that there is a good deal of new physics beyond the standard model. 
In cosmology, the SM does not explain dark matter, inflation and dark energy. On the side 
of particle physics, the recent discovery of flavor conversion of solar, atmospheric, reactor, 
and accelerator neutrinos has established conclusively that neutrinos have nonzero mass and 
they mix among themselves much like quarks, thereby providing the first evidence of new 
physics beyond the standard model. This discovery has triggered enormous theoretical ac- 
tivity attempting to uncover the nature of this new physics. These attempts include further 
developments of the already existing mechanisms and theories such as the grand unified 
theories (GUT's) and the appearance of new ideas and approaches. Various aspects of these 
developments have to some extent been covered in several recent reviews jlj]. The present 
review is an update that focuses on what we think are the most important perspectives on 
the most plausible and natural scenarios of physics beyond the standard model. 

The main points stressed in this review are as follows: 

1) . After all recent developments, the seesaw mechanism with a large scale of B — 
L (baryon minus lepton number) violations still appears to be the most appealing and 
natural mechanism of neutrino mass generation. However, it is not excluded that some 
more complicated version of this mechanism is realized. 

2) . Grand unification (plus supersymmetry in some form) still appears to be the most 
preferable (plausible) scenario of physics that naturally embeds the seesaw. 

3) . The seesaw - GUT scenario does not provide a complete understanding of the masses 
and mixing of neutrinos, nor the masses and mixing of other fermions (in other words, the 
flavor structure of the mass matrices). Some new physics in addition to this scenario seems 
essential. In this connection two issues are of great importance: 

- possible existence of new symmetries that show up mainly (only?) in the lepton sector; 
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and 

- understanding the relation between quarks and leptons. 

To uncover this new additional element (s) of theory, the bottom- up approach is extremely 
important. 

4). Alternative mechanisms and pictures are not excluded, although they have not 
reached the same level of sophistication as the GUT approach and appear less plausible. 
Nonetheless, it is quite possible that they play a role in the complete picture, e.g., as the 
source of new neutrino states, as some element in addition to the "seesaw - GUT" scenario, 
or as some element of physics at a level deeper than or above the GUT scale. Flavor structure 
can appear as a result of compactification of extra dimensions. 

This review presents a detailed discussion of these statements and current arguments 
in favor of this possibility. In sect. 2 we introduce the main notions on neutrino mass 
and mixing and summarize available experimental results. We proceed with the bottom-up 
approach in the following three sections. In sect 3. we analyze results on neutrino mixings 
and the mass matrix. We address the question of particular neutrino symmetries in sect 4., 
and in sect. 5 we consider possible relations between quarks and leptons. In the rest of the 
review the top-down approach is presented. We study properties of the seesaw mechanism 
in sect 6. Neutrino masses in the grand unified theories are discussed in sect. 7. We consider 
achievements and limitations of this picture. The need for flavor symmetries in addition to 
GUTs is described. We also discuss alternative mechanisms of neutrino mass generation and 
the alternative scenarios for the big picture. In sec. 8, we discuss sterile neutrinos, their 
imphcations and possible theories. Sec. 9 contains the conclusion. 

II. MASSES, FLAVORS AND MIXING 
A. Neutrinos and the Standard model 

According to the standard model (SM), left-handed neutrinos form the EW doublets L 
with charged leptons, have zero electric charge and no color. The right-handed components, 
vr, are not included by choice. 

The masslessness of the neutrinos at the tree level in this model owes its origin to the 

fact that there are no right-handed neutrinos. This result holds not only to all orders in 
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perturbation theory but also when nonperturbative effects are taken into account, owing 
to the existence of an exact B-L symmetry in the model even though B+L is violated by 
weak sphaleron configurations. It would therefore appear that nonzero neutrino masses must 
somehow be connected to the existence of right-handed neutrinos and/or to the breaking of 
B-L symmetry both of which imply new physics beyond the standard model. 

Vanishing conserved charges (color and electric) distinguish neutrinos from other fermions 
of the standard model. This difference leads to several new possibilities for neutrino masses 
all of which involve new physics: 

(i) The neutrino masses could be the Majorana type thereby breaking L by two units. 

(ii) Neutrinos have the possibility to mix with singlets of the SM symmetry group, in par- 
ticular, singlet fermions in extra dimensions. 

The main question is whether these features are enough to explain all the salient 
properties of the neutrino masses and mixing observed in experiment. 



One way that the neutrino mass can be generated even if the SM particles are the only 
light degrees of freedom requires that the condition of explicit renormalizability of the theory 
should be abandoned. Indeed, the non-renormalizable operator [2I 

^(L.HfiL.H), 2,j = e,/i,r, (1) 

where H is the Higgs doublet Xij are the dimensionless couplings and M is the cutoff scale, 
after the EW symmetry breaking generates the Majorana neutrino masses 

(^) 

The operator breaks L and B — L quantum numbers. One may think that this operator is 
generated by some gravitational Planck scale effects, so that M ~ Mpi and Xij ~ 1 [3]. In 
this case however, rriij ~ 10^^ eV are too small to explain the observed masses (although 



such a contribution can still produce some subleading features Therefore, new scales of 
physics below Mpi must exist to give the desired mass to neutrinos. The operator in eq.([l]) 
can appear after integrating out some new heavy degrees of freedom with masses M <^ Mpi. 

Another important conclusion from this consideration is that the neutrinos can get rel- 
evant contributions to masses from all possible energy/mass scales M from ~ 1 eV to the 
Planck scale. If there are two or more differentsubstantial contributions (from different scales 



and different physics) to the mass, the interpretation of results can be extremely difficult. 
In this case one can argue that the "coincidence" of sizes of contributions looks unnatural. 
However, this will be not the only place where we face the coincidence problem. 

B. Right-handed neutrinos, neutrino mass and seesaw. 

Let us consider possible extensions of the standard model which can lead to non-zero 
neutrino masses. 

1) . If right-handed neutrinos exist (we consider the conceptual implications of their 
existence in sect. 6.2) one can introduce the Yukawa coupling 

Y^LHur + h.c. (3) 

which after the electroweak (EW) symmetry breaking leads to the Dirac neutrino mass 

mD = Y,{H). (4) 

The observed neutrino masses would require < 10^^^ — 10^^^. If ur is the same type of 
field as right-handed components of other fermions, such smallness looks rather unnatural. 
However, if the Dirac mass is formed by coupling with some new singlet fermion, S, beyond 
the usual fermion family structure, possible new symmetries associated to S and/or z/^ can 
suppress Y^. In this case Y^, appears as the effective coupling: Yi^ ~ (f^/M)", where f 5 ^ M 
are the scales of some new interactions and new symmetry breaking, and n is some integer 
determined by the charges of the fields. 

2) . The right handed neutrinos are allowed to have Majorana masses 

MrvIC-^vr + h.c., (5) 

where C is the Dirac charge-conjugation matrix. Because ur are singlets under the SM 
gauge group the Mr can appear as a bare mass term in the Lagrangian or be generated by 
interactions with singlet scalar field a (so that Mr — > f„a in eq.Q): 

MR = Ua), (6) 

where (a) is the vacuum expectation value (VEV) of a. The latter possibility is realized if, 
e.g., vr is a component of a multiplet of an extended gauge group. 
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3). The left handed neutrinos can also acquire the Majorana masses m^. The corre- 
sponding mass terms have the weak isospin 1=1 and violate lepton number by two units. 
Thus, they can be generated either via the non-renormalizable operators in eq.([T]) with two 
Higgs doublets and/or by coupling with the Higgs triplet A: 

fAL'^LA + h.c. (7) 

The non-zero VEV of A then gives rriL = /a (A). 

In the case of three neutrino species the mass parameters rriD, Mr and rriL should be 
considered as 3 x 3 (non-diagonal) matrices. In general all these mass terms are present. 
By introducing the charge-conjugate left handed component A^^^, = {vr)^ the general mass 
matrix in the basis {vl-, Nl), can be written as 



( T 



rriD Mr 



(8) 



The eigenstates of this matrix are the Majorana neutrinos with different Majorana masses. 
The Dirac mass term mixes the active neutrinos with the sterile singlet states N^. 
The matrix has several important limits. 

Suppose niL = 0. The Yukawa couplings eq.Q are expected to be of the same order 
as the charged-fermion couplings. Because the Nls are singlets under the SM gauge group, 
their Majorana masses unlike the masses of the charged fermions, are not constrained by 
the gauge symmetry and can therefore be arbitrarily large, i.e. Mr rrijj. In this case 
the diagonalization of the mass matrix in eq.([H]), leads to an approximate form for the mass 
matrix for the light neutrinos m^: 

M, = -mlM^'niD. (9) 

As noted, because Mr can be much larger than mD one finds that rrii, -C me,u,d very 
naturally, as is observed. This is known as the seesaw (type I) mechanism [sl and it provides 
a natural explanation of why neutrino masses are small. 

If elements of the matrix rriL are non-zero but much smaller than the other elements of 
Mjy, we can write the resulting light neutrino mass matrix in the form 

Mu = rriL - m^M^i^mD. (10) 
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We refer to this as mixed seesaw [6|, [7|] and when the first term dominates, we refer to it as 
a type II seesaw. 

In the matrix eq.([S]), elements of both and Mr may have magnitudes much smaller 
than those of m/j. In this case, the neutrinos are predominantly Dirac type, with a small 
admixture of Majorana mass. We refer to this case the pseudo-Dirac [8]. 

The main question we discuss subsequently is whether the see-saw mechanism can explain 
all the observed features of neutrino mass and mixing. 



C. Flavors and mixing 

The electron, muon and tau neutrinos, z/g, v^, - particles produced in association 
with definite charged leptons: electron, muon and tau correspondingly are called the flavor 
neutrinos. For example, the neutrinos emitted in weak processes such as the beta decay or 
pion decay with electron or muon are termed the electron or muon neutrinos. In the detection 
process, it is the flavor neutrino that is picked out because the detectors are sensitive to the 
charged lepton flavors such as (e, /i,r). In the SM neutrinos v^.^ i/^, and are described by 
the flelds that form the weak doublets (or weak charged currents) with charged lepton flelds 
of deflnite mass: 

j^ = lY{X-^^)vi, / = e,/i,r. (11) 

Note that phenomenologically deflned flavor states as states produced in certain weak 
processes, may not coincide precisely with theoretical flavor states - the states produced by 
the flelds from certain weak doublets. The difference can appear because of neutrino mixing 
with heavy neutral leptons that cannot be produced in low-energy weak processes owing to 
kinematics. In fact, this situation is realized in the seesaw mechanism. However, in standard 
situation, the admixture and difference of the states is negligible. 

Flavor mixing means that the flavor neutrinos z/q, [a = e, fi, r) do not coincide with 
neutrinos of deflnite mass z/j [i = 1,2,3). The electron, muon and tau neutrino states have 
no deflnite masses but turn out to be coherent combinations of the mass states. The weak 
charged current processes mix neutrino mass states. 

The relationship between the flavor flelds, uj = (z/g, z/^, Ur) and the mass flelds u"^ = 
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(z/i, z/2, z/3) can be written as 

= UpMNS^, (12) 

where Upmns is 3 x 3 unitary matrix termed the Pontecorvo - Maki - Nakagawa - Sakata 

n □ 

lepton mixing matrix [9|, [lOl]. The fields z/J = (z/g, t'^, z/,-) form the fiavor basis. Inserting 
eq. (fT2!) into eg. (Ill I) we can write the weak charged currents as 

J^^ = lYil-^,)UpMNSl^- (13) 

Thus, the lepton mixing matrix connects the neutrino mass fields and charged lepton fields 
with definite mass in the weak charged currents. Neutrino mass states are the eigenstates 
of the Hamiltonian in a vacuum, and we refer to the mixing in eq. fll2p as vacuum mixing. 
Vacuum mixing is generated by the non-diagonal mass matrices. 

In addition to the fiavor Uf and mass bases let us introduce the "symmetry" basis (z/, /) 
- the basis in which the underlying theory of the fermion masses is presumably formulated. 
This can be some fiavor symmetry, or GUT, or some dynamical principle, or selection rule 
originating from string theory. We do not know this basis a priori, and in fact, its identifi- 
cation is one of the key problems in the bottom-up approach. 

In the symmetry basis, both the charged lepton and the neutrino mass matrices are in 
general nondiagonal (although models exist in which the symmetry basis coincides with the 
fiavor basis). The mass terms of the Lagrangian can then be written as 

Cm = ulC-'M^UL+'hMiln + h.c. (14) 

We have assumed that neutrinos are Majorana particles. (Notice that in the fiavor basis 
(z//, /) the mass matrix of charged leptons is diagonal, so the existence of mixing implies 
that the mass matrix of neutrinos should be non-diagonal in this basis.) We diagonalize the 
matrices in eq. (ll4l) as 

UjMM. = Mt, UiMiVi^ = Mf, (15) 

where = diag{mi, 1712,1713) and = diag{me,m^,mT-) are the diagonal matrices and 
the rotation matrices, Uu, Ue, Ve, connect the symmetry fields with the mass fields: u = U^u, 
II = UiIl, Ir = ViIr. Plugging these relations into the charged current we obtain 

= [^^.(l - 75)z> = /Y(l - 75)[//[/,z.. (16) 
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Thus, the physical neutrino mixing matrix is given by 

UPMNS = U}U,. (17) 
It is convenient to parameterize the mixing matrix as 

UpMNS = f/23(^^23)f/l3(^13, 5)f/l2(^12)/<^, (18) 

where the Uij are matrices of rotations in the ij plane by angle 6'jj, 5 is the Dirac CP- 
violating phase attached to 1-3 rotation. In the case of Majorana neutrinos sometimes the 
mixing matrix is defined as f/pMA^s = UpmnsI^), where = diag{l, 6**^^ e**^^) is the diagonal 
matrix of the Majorana CP-violating phases. The latter can be incorporated into the mass 
eigenvalues which can then be considered as complex parameters. 

Note that in the first approximation the above parameterization allows us to connect im- 
mediately the rotation angles with physical observables 6*23 = Oatm, is the angle measured in 
the atmospheric neutrino oscillations; 612 = dsoi is the angle determined from solar neutrino 
studies, and ^^13 = Ochooz is the angle restricted by the reactor experiment CHOOZ. 



D. Experimental results and global fits 

We use the results of global analysis of the neutrino data published through the middle of 
2006. The analysis assumed that (i) there are only three mixed, active neutrinos; (ii) CPT 
is conserved, so that masses and mixing angles in the neutrino and antineutrino channel 
coincide; and (iii) neutrino masses and mixings have pure "vacuum origin": that is, due 
to the interaction with Higgs field(s) that develop a VEV at a scale much larger than 
neutrino mass. We comment below on possible changes when some of these assumptions are 
abandoned. 

The parameter space includes the oscillation parameters: mass-squared differences 
Amfj = — rrij, mixing angles % as well as the Dirac CP-violating phase S. Non- 
oscillation parameters are the absolute mass scale which can be identified with the mass of 
the heaviest neutrino and two Majorana CP-violating phases. 

The experimental results used in the analysis can be split into three sectors: 



1). Solar neutrinos ll| and the reactor experiment KamLAND 121] are sensitive to 



mainly Am^i and 612 (solar sector). The 1-3 mixing, if not zero, may give subleading effect. 
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2) . Atmospheric neutrino studies [13|, K2K [Q and MINOS [? ] accelerator experiments 
are sensitive to Am^^ and 6*23 (atmospheric sector). The solar parameters Am^i and 6*12 as 
well as 6*13 give small subleading effects. 

3) . The CHOOZ experiment [l5| gives a bound ^13 as a function of Aml^. 
The physical effects involved in the interpretation are 

- vacuum oscillations [o], lo|, IgI (atmospheric neutrinos - main mode, K2K, MINOS and 
CHOOZ); 

n 

- the Mikheyev-Smirnov-Wolfenstein (MSW) effect - the adiabatic conversion |17l . 
(conversion of solar neutrinos in the matter of the Sun; at low energies solar neutrinos 
undergo the averaged vacuum oscillation with small matter effect); and 

n n 

- oscillations in matter |17l . Il8l | (oscillations of solar and atmospheric neutrinos in the 
matter of Earth). These oscillations produce sub- leading effects and have not yet been 
estabhshed at the statistically significant level. 



Let us consider the results of global analysis from [19|, |20[ |. 

The mass splitting responsible for the dominant mode of the atmospheric neutrino oscil- 
lations equals 



|Am22| = (2.4 ± 0.3) ■ lO'^ eV^, (la). 



(19) 



The sign of the mass split determines the type of mass hierarchy: normal, > 0) 

inverted, Aml2 < 0, and it is not identified yet. The result of eq. (fT9|) allows one to get a 
lower bound on the heaviest neutrino mass: 



ruh > Vl^^isl > 0-04 eV, (2a), 



(20) 



where m/j = for the normal mass hierarchy, and m/j = mi ^ m2 for the inverted hierarchy. 

A much smaller mass-squared splitting drives the solar neutrino conversion and oscilla- 
tions detected by KamLAND 



Am22 = (7.9 ± 0.4) ■ 10"^ eV^ (la). 



(21) 



In the case of the hierarchical mass spectrum that would correspond to m2 ~ 0.009 eV. The 
ratio of the solar and atmospheric neutrino mass scales. 



Am|i 
Amii 



0.033 ± 0.004, 



(22) 
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FIG. 1: The best-fit values and the allowed regions of lepton mixing angle 612 at different confidence 



levels determined by different groups. From SNO (from ref. 



ui), sv 



19| and Bari [20] . Shown are 



predictions from quark-lepton complementarity (QLC) and tri-bimaximal mixing (see sees. 4,5). 
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FIG. 2: The best-fit values and the allowed regions ofsin^ ^23 at different confidence levels deter- 
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mined by different groups: SK 
5) and sensitivity limit of T2K experiment 



m, Bari 



201 1 . Shown are expectations from QLC (sec. 



gives a lower bound on the normal mass hierarchy 

— > v^ = 0.18 ±0.01. 



(23) 



In figs. [U [3] we summarize results of determination of the mixing angles obtained by 
different groups. We also show some theoretical benchmarks that will be discussed below. 
The best-fit value of the 1-2 mixing angle and la error equal 



= 33.9° ± 1.6°, 



(24) 
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FIG. 3: The best-fit values and allowed regions of sin^ ^13 at different confidence levels determined 
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23] and T2K 21|]. 



or sin^ 612 



U.OiO_Q 



[see fig. [T]). It deviates from maximal mixing, sin^ ^12 = 0.5, by 



approximately (6 — 7)cr. 

The 2-3 mixing is in agreement with maximal mixing, 623 = it /A (fig. [2]). A shift of the 
best-fit point from 7r/4 to smaller angle s appea rs when the effect of the 1-2 mass splitting 



and mixing is included in the analysis 



20 



23j. According to |23|, sin^^ss = 0.47, and a 



slightly larger shift, sin^ ^23 = 0.44, follows from analysis 20]. Thus, the deviation from 



maximal mixing can be quantified as 

D23 = 0.5 - sin^ 9. 



23 



0.03 - 0.06. 



(25) 



The shift is related to an excess of the so-called e-like atmospheric neutrino events in the sub- 
GeV range. The excess can be explained by the oscillations driven by the "solar" oscillation 
parameters and it is proportional to the deviation D23 2J]. The experimental errors (the 
lower bound sin^ 26'23 > 0.92) still allow substantial deviations from maximal mixing: 

D23/sin2^23~0.4 {2a). (26) 

Results on the 1-3 mixing are consistent with zero ^13 (fig. [3]). A small, nonzero best-fit 
value of sin^ 6*13 from the analysis [20|] is related to the angular distribution of the multi-GeV 
e-like events measured by Super-Kamiokande 13]. The most conservative 3a bound is [2^ 

sin^ ^13 < 0.048, (3a). (27) 

In the first approximation, the pattern of lepton mixing has been established. There are 
two large mixings: The 2-3 mixing is consistent with maximal mixing, 1-2 mixing is large 
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but not maximal, and 1-3 is small and consistent with zero. Further precise measurements 
of the mixing angles and, in particular, searches for the deviations of 1-3 mixing from zero 
and of 2-3 mixing - from maximal, are crucial to understanding the underlying physics. 
The figures show the accuracy required to make important theoretical conclusions and the 
potential of next generation experiments. 



Information on non-oscillation parameters has been obtained from the direct kinemati- 
cal measurements, neutrinoless double beta decay and cosmology. The effective Majorana 
mass of the electron neutrino - the ee-element of the neutrino Majorana mass matrix, mee, 
determines the life time of the neutrinoless double beta decay: T(2/50i/) oc m~^. In terms of 
masses and mixing parameters, it can be written as 



nif. 



k 



i(t>{k) I 



k 



i<t>{k) 



(28) 



where (/)(fc) is the phase of the k eigenvalue, Uek its admixture in z/g, and - the lightest 
neutrino mass. Fig. H] from |l^ summarizes the present knowledge of the absolute mass 
scale. The regions allowed by oscillation results in the plane of and mi - the mass of 
the lightest neutrino probed by the direct kinematical methods and cosmology are shown. 
The two bends correspond to the normal and inverted mass hierarchies. For a given mi, the 
range of rriee is determined by variations of the Majorana phases and uncertainties in the 
oscillation parameters. 

The Heidelberg-Moscow experiment gives the best present bound on mee'- f^iee < (0.35 — 
0.50) eV 251]. Part of the collaboration claims an evidence for a positive signal 26j which 
would correspond to niee ~ 0.4 eV. If this positive signal is confirmed and if it is due to the 
exchange of light Majorana neutrinos, the neutrino mass spectrum is strongly degenerate: 
mi ^ m2 ~ m3 = mo 27 1. 

The cosmological observations put a bound on the sum of neutrino masses 

< 0-42 eV (95% C.L.) |28j (see also [291) which corresponds to mo < 0.13 eV in 
ihe case of a degenerate spectrum. An even stronger bound, J2i=i''T^i < 0.17 eV (95% C.L.) 
3o| was established after publication of WMAP3 results. This limit disfavors a strongly 
degenerate mass spectrum. Combining the cosmological and oscillation Eq. fl20p bounds, 
we conclude that at least one neutrino mass should be in the interval 



m ~ (0.04 - 0.10) eV (95% C.L.). 



(29) 
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FIG. 4: The 99% C.L. range for rriee as a function of the hghtest neutrino mass for the normal 
(Am23 ^ 0) ^^^^ inverted (Am23 < 0) mass hierarchies. The darker regions show the allowed range 
for the present best-fit values of the parameters with negligible errors; from [2]. 



Direct kinematic measurements give the weaker limit m < (2.0 — 2.2) eV 3l|]- The 



planned experiment KATRIN 32| is expected to improve this limit down to ~ 0.2 eV. 

How robust are these results? Can we expect some substantial change in this picture in 
future? There are three types of effects (in fact, related to lifting of assumptions made in 
the analysis) that can influence the interpretation of the present neutrino results: 

1) . Possible existence of new neutrino states - sterile neutrinos. If these states are light, 
they can directly (dynamically) influence the observed effects used to determine neutrino 
parameters. 

2) . Possible presence of non-standard neutrino interactions can change values of the 
extracted neutrino parameters. 



3). Interactions with light scalar fields 33| that can produce the soft neutrino masses 



that depend on properties of the medium. These masses may also change with time and be 
related to dark energy in the universe j34|. 

At present, however there are no well-established results that indicate deviation from the 
"standard" 3z/ mixing and standard matter interactions. 
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NORMAL INVERTED 



FIG. 5: Neutrino mass and flavor spectra for tlie normal (left) and inverted (right) mass hierarchies. 
The distribution of flavors (colored parts of boxes) in the mass eigenstates corresponds to the best- 
fit values of mixing parameters and sin^ ^13 = 0.05. 

There are various tests of validity of the standard picture and of the theory of neutrino 
conversion. Different sets of data confirm each other. Consistent interpretation of whole 
bulk of various results in terms of the vacuum masses and mixing provides further 
confidence. The fit of the data is not improved substantially with inclusion of new states 
and non-standard interactions and, if exist, they may produce sub- leading effects only. One 
can perform, e.g., a global fit of neutrino data considering the normalization of the matter 
potential as free parameter. According to [20] the best-fit value of the potential is close to 
the standard one. In this way one tests not only the validity of the refraction theory for 
neutrinos, but also the consistency of the whole picture. 

E. Mass and flavor spectrum 

Information obtained from the oscillation experiments allows us to partially reconstruct 
the neutrino mass and fiavor spectrum (Fig. [5]). There are four unknowns: (i) the admixture 
of i>e in 1/3 described by t/es! (ii) the type of mass spectrum - hierarchical; non-hierarchical 
with a certain ordering; or degenerate, which is related to the value of the absolute mass 
scale, mi; (iii) the type of mass hierarchy (ordering): normal, or inverted; and (iv) the 
CP-violating phase 6. Determining these unknowns is the goal of future phenomenological 
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and experimental studies. 

There are some weak indications in favor of a normal mass hierarchy from supernova 
SN1987A data. However in view of small statistics and uncertainties in the original fluxes 
it is not possible to make a firm statement. 

As is clear from the fig. IH future high-sensitivity measurements of the effective mass rriee 
can allow one to establish the hierarchy. The bound rriee < 0.01 eV will exclude the inverted 
mass hierarchy and also degenerate mass spectrum if neutrinos are Majorana particles. 
Future detection of a galactic supernova can also contribute to the determination of the 
type of mass hierarchy and 1-3 mixing 35 1. 



F. Toward the underlying physics 

What is behind all these observations? To uncover the underlying physics two general 
strategies are invoked: 

1. The Bottom- Up approach is essentially an attempt to uncover the underlying physics 
starting from observations. The strategy is to (i) reconstruct the neutrino mass matrix in 
the flavor basis using the information available on masses and mixings {Amjj, 6ij, irtee)] (ii) 
take into account the RGE effect and obtain the mass matrix at the scale of new physics, 

(iii) search for the symmetry basis in which flavor or some other symmetry is realized; and 

(iv) identify the symmetry and mechanism of symmetry violation, if needed, as well as the 
underlying dynamics. This is our standard way of understanding things, but it does not 
preclude that the explanation will require something new. 

2. Top-Down approach starts with a general unifled theory framework, be it a GUT, TeV 
scale theory or extra dimension theory which has motivation outside neutrino physics - and 
uses it to make predictions for neutrinos, i.e., to go from the big picture to the observed 
properties of neutrinos. 

At some point these two approaches should merge. Both approaches are needed. It seems 
difficult to uncover the underlying picture by just moving from observations, and some a 
priori ideas (or context) are needed to relate neutrinos with other physics. This opens the 
possibility of using results from other areas {e.g., collider experiments). Thus, we need the 
big picture. However, working solely within the big picture, one can miss some important 
elements, which is where the bottom-up approach can help. The results of the bottom-up 
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approach are considered in sect. 3-5. The top-down approach is developed in sect. 6-7. 



III. MIXING AND MASS MATRIX OF NEUTRINOS 

Analyzing results on neutrino mixing and mass matres one can (i) search for some par- 
ticular features in the data, such as empirical relations, equalities, hierarchies of elements, 
zeros, etc.; (ii) identify possible dominant structures in the mixing and mass matrices (the 
idea being that matrices can have a structure as "lowest order plus small corrections" which 
in turn can correspond to some dominant mechanism plus sub-dominant effects); (iii) search 
for simple parameterizations in terms of small number of parameters; and (iv) present ma- 
trices in terms of powers of some small quantity, etc.. All these may hint at the underlying 
physics. 



A. Properties of neutrino mixing matrix 



Let us first to analyze the mixing matrix. Maximal 2-3 mixing, large 1-2 mixing and 
small 1-3 mixing indicate that the following matrices can play some important role {e.g., be 
the lowest-order mixing matrices). 



1). The bimaximal mixing matrix 36|: 



TT — TTinTj-m 
'^bm — '-^23'-^ 12) 



(30) 



where and are the matrices of the maximal (vr/4) rotations in the 2-3 and 1-2 
subspaces respectively. Explicitly, we have 



U, 



bm 



( V2 V2 \ 



V 



(31) 



-1 1 V2 
1 -1 72/ 

The identification Upmns = Ubm is not possible owing to strong deviation of the 1-2 mixing 
from maximal. However, Ubm can play the role of a dominant structure. In the latter 
case, the correction can originate from the charged lepton sector (mass matrix), so that 
UpMNS = U'Ubm- Suppose U' ~ Ui2{(y), where a ~ 0{9c) in analogy to quark mixing. 
Then U' simultaneously generates deviations of the 1-2 and 2-3 mixing from maximal as 
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well as non-zero 1-3 mixing which are related as 



sin 



13 



sm a 
~V2 



and 



sin &■ 



23 - (l - sin^ 6*13 



or D23 ~ sin^ ^13. Confirmation of this equality will be very suggestive. 
2). The tri-bimaximal mixing matrix 37 1 



t/tfem = f/2^f/i2(^i2), sin2^i2 = l/3, 



(32) 



or explicitly, 

/ 2 72 ^ 

Utbm = ^ -1 V2 V3 (33) 
\ 1 -v^ V3y 

is in good agreement with data including 1-2 mixing. Here V2 is trimaximally mixed: In the 
second column the three flavors mix maximally, whereas is bi-maximally mixed. Mixing 
parameters turn out to be simple numbers such as 0, 1/3, and 1/2 and can appear as the 
Clebsch- Gordon coefficients. 

The bi-maximal and tri-bimaximal matrices can be considered as matrices in the lowest 



order of some approximation. Then one can introduce, 
of the true matrix from these lowest-order structures 



parameters that describe deviations 
The matrices Uhm and Utbm reveal 



certain symmetries, and the deviation parameters may describe effects of the violation of 
these symmetries. 



B. Reconstructing neutrino mass matrix 

The mass matrix is probably more fundamental than mass eigenvalues and mixing angles 
because it combines information about masses and mixings. Dynamics and symmetries can 
be realized in terms of mass matrices and not their eigenstates and eigenvalues. However, it 
is possible (and models of this type exist, see below) that symmetry immediately determines 
the mixings but not the masses which are left as free parameters. 

As mentioned above, the first step in the bottom-up approach is the reconstruction of the 
mass matrix in flavor basis. Note that in the case of Majorana neutrinos, the elements of this 
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mass matrix are physical parameters. They can be measured directly, e.g., in neutrinoless 
double beta decay and, in principle, in other similar processes. 

The answer to the question what is more fundamental: mass matrices or observables 
(Am^, 6), may depend on the type of mass spectrum. In the case of a hierarchical spectrum, 
the observables are imprinted visibly into the structure of the mass matrix. In contrast, 
for the quasi-degenerate spectrum they are just very small perturbations of the dominant 
structure determined by the non-oscillatory parameters: the absolute mass scale and the 
Majorana CP-violating phases. The oscillation parameters can originate from some small, 
in particular radiative, corrections. 

In the flavor basis the mass matrix of the charged leptons is diagonal and therefore the 
neutrino mass matrix is diagonalized by f/j, = Upmns- Consequently, according to eq. (fTH|) 
the neutrino mass matrix in the flavor basis can be written as 

= U}.,,^sMtU^pMNS. (34) 

where 

Mt = diagimi, m^e-'^''^^ mge-^*'^^). (35) 

Here ^are the Majorana phases and we take here 0i = 0. Results of reconstruction 



show [39[ that a large variety of different structures of mass matrices is possible, depending 
strongly on the unknown mi, type of mass hierarchy and Majorana phases. The dependence 
on sin and 5 is relatively weak. This means huge degeneracy of mass matrices now 
and perhaps even in the far future because, in reality, it is not possible to measure all the 
parameters including CP-violating phases. Variations of one Majorana phase (even if all 
other parameters are known) can lead to a strong change in the structure. Nature would 
have to be very "collaborative" for us to determine the mass matrix. Or we may uncover 
some principle that allows us to predict the mass matrix, which we can then check by certain 
precision measurements. 

C. Extreme cases 

To give some idea about various possibilities, we present simple parameterizations of the 
neutrino mass matrix in the flavor basis for three extreme cases: normal mass hierarchy, 
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inverted mass hierarchy and degenerate spectrum. 



1). Normal mass hierarchy (mi <^ m2 -C ^3). The mass matrix indicated by data can 
be parameterized as 

/ de be ae \ 



31 



(36) 



be 1 + ce -1 
\ae -1 l + e/ 

where a, b, c, d are complex parameters of the order one, and e is essentially the ratio of the 
solar and the atmospheric mass hierarchies: e = 2^yrXF{a, c, b, d), and F(a, c, b, d) ~ 0(1). 
A salient feature of this matrix is the dominant — r block. With the present accuracy 
of measurements of the parameters, a sharp difference between dominant and subdominant 
elements may not exist and some moderate hierarchy of elements with the unique expansion 
parameter 0.6 - 0.7 is realized 391 ]. 

An important property of the above mass matrix is that in the limit of a = 6 and c = 1, it 
is symmetric under /i — r interchange and one gets maximal 2-3 mixing and zero 1-3 mixing. 



I.e. a 



23 



? and 9 



13 







40l |. This symmetry should of course be approximate one, as the 



masses of muon and tau leptons are different. Any resulting /i — r breaking must therefore 



reflect itself in nonzero 61^ and D23 with the two connected to each other [4l|. For example, 
if a = 6 and c 7^ 1, i.e., the symmetry is broken in the dominant block, the induced 6'i3 and 
D23 are given by 



sin 9 



13 



-be\c-l)/A, D23 = e(c-l)/4, 



and therefore are related as 



tan6'i3 = beD 



23- 



(37) 



(38) 



Furthermore, ^13 ^ 1^23 • In contrast, if the symmetry is broken in the sub-dominant block 
only: a ^b but c = 1, the situation is opposite: ^13 ^ D23, ^-c 



sin^i 



a — b 
~V2 



D 



23 



-e . 



(39) 



Thus, measurements of ^13 and 1^23 will provide an important probe of the mass matrix 
structure. Note that when a = b = d and c = 1, we get the tri-bi-maximal mixing pattern. 
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2). Inverted mass hierarchy (mi ^ m2 ^ ^713). The structure of the mass matrix in this 
case depends strongly on the Majorana CP-violation phase. An approximate form of the 
mass matrix in the case of opposite CP-parities of ui and z/2 is 



where e -C 1. In the limit e — > 0, this mass matrix has the symmetry — L^^ — |I42 |. 
In the symmetry limit one has Am^2 = 612 = 7r/4. Furthermore, if an additional 
/i — r exchange symmetry is imposed on this mass matrix, the atmospheric mixing angle 
also becomes maximal. 

The breaking of — — symmetry leads to nonzero Am^2 ^'^d deviation of 612 
from maximality. It has, however, been difficult in general, although not impossible, to 
accommodate the observed "large" departure from maximality of 612 using "small" breakings 
of Lf. — — Lt- symmetry. One needs as much as 40% breaking to fit data. 

In the case of the same CP-parities of mass eigenstates the mass matrix has completely 
different form with interchange of dominant - sub-dominant elements in eq.( HOl) . This 
illustrates strong dependence of the matrix structure on the unknown CP violating phases. 

3). Degenerate spectrum: mi ^ m2 ~ = ttlq. 

Here also the structure of mass matrix depends strongly on the CP-violating phases. Two 
possibihties are particularly interesting: 

(i) If the relative phases between the mass eigenvalues are zero (27r/c), the mass matrix 
is close to the unit matrix 



ze c 





c ye de 
y s de xe J 



(40) 



M 



V 



mo/ + 5M, 



(41) 



where 8M <^ mo is the matrix of small corrections, 
(ii) In the case of opposite CP-parities of z/2 and v-i, 



/l o\ 



M 



1 



+ 5M. 



(42) 



Vo 1 oy 
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Both matrices give rriee = ttlq and can explain the Heidelberg-Moscow positive result. 
Theoretical understanding of such a situation would require underlying symmetries that 
guarantee mass degeneracy. The basic strategy here is to consider symmetries that have a 
three-dimensional representation to which the three lepton doublets of the standard model 
can be assigned and then design a Higgs sector that will lead to a quasi-degenerate neutrino 



spectrum. A list of such symmetries includes [43(], 5*0(3) [44], and A4, j45|. 

In this connection an interesting possibility is a theory with mixed seesaw, where the 
type II contribution gives the dominant quasi-degenerate term itiqI or the first term in 
eq. ( l42l) with the conventional type I contribution giving mass splittings and mixings: 5M = 
—m/^M^^rriD. A very generic way to see how these models could explain observations is 
as follows: In a quark-lepton unified theory such as an SO(IO) model, we would expect the 
Dirac mass term for the neutrinos to have a hierarchical pattern for its eigenvalues so that 
roughly speaking, the atmospheric and solar mass differences will be given by 

respectively, roughly similar to observations. 

The matrices of eqs. (14111421) offer various possibilities that relate the degeneracy of the 
spectrum with large or maximal mixing. The matrix of eq. (H2l) leads immediately to maximal 
2-3 mixing. Nearly maximal mixings are generated by small off-diagonal elements in eq.( HTi) . 



D. Mass matrices with texture zeros 



Another approach in analyzing possible mass matrices is to see if some elements can be 
exactly zero or equal each other. This may also uncover dominant structures and certain 
underlying symmetries. This approach allows one to reduce the number of free parameters 
and therefore can lead to certain predictions. Recall that the Majorana mass matrix for 
three neutrinos has six independent elements. 

Mass matrices with different numbers of zeros and with zeros in various places of matrix 
have been considered. Two of the cases discussed widely in the literature are textures with 
three zeroes and two zeros. It is easy to be convince that the three zeros cannot be along 
the off-diagonal entries nor be in any of the 2x2 submatrices and still give a fit to already 
known data. In the former case all mixings vanish and in the latter case, one cannot satisfy 
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the requirement from observation that Am^2 ^ ^'^as if ^23 and 612 are large, as observed. 
The case when all zeros are along the diagonal [46] (or two along the diagonal and the 
third is off diagonal) is more subtle because now one can satisfy the requirements of large 
solar and atmospheric mixings as well as Am^g ^ ^f^ls- However, in this case, there are 
only three (real) parameters which can be determined from Am^2 Am23 ^'^'^ ^23- One then 
predicts a value for the solar mixing angle, sin^ 2^12 = 1 — rA/16, which is incompatible with 
observations. 

As far as textures with two zeros are concerned, they have five free parameters: four real 
parameters and a complex phase and are therefore interesting candidates for neutrino mass 
matrices {4^]. These have been analyzed to give their characteristic predictions. There are 
seven different (out of fifteen) possibilities that are currently in accord with data and make 
predictions for various parameters such as neutrinoless double beta decay and ^^13. As an 
interesting example, consider the matrix 

/ X\ 



(43) 



X X 
\X X X J 

where X indicate non-zero entries. This leads to a hierarchical mass matrix with the pre 
diction 



sin^^ 



13 



tan^ 6*1 



/12 - cot^ 9i2 

and zero amplitude for neutrinoless double beta decay. 
Another possible texture is 

fx X 0\ 



0.01 



(44) 



(45) 



X X 

V XX/ 

This produces a degenerate mass spectrum with an effective mass in neutrinoless double 
beta decay exceeding 0.1 eV. 

Such an approach should be taken with some caution: (i) for instance, it is not clear why 
the zeros appear in the flavor basis? (ii) there is the possibility that interesting symmetries 
may not correspond to zeros in the flavor basis, and (iii) finally, the zeros may not be exactly 
zeros, in which case we are unlikely to learn much from these exercises. In any case, such 
situations should be considered on equal footing with cases in which various elements of the 
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neutrino mass matrices are related or simply equal, as in the case of /i — r symmetric models 
discussed in sec. 4. 



E. Anarchy approach 

Although the limiting examples of mass matrices described above may contain interesting 
hints of symmetries, it is quite possible that one is far from these cases. In a large part of 
the parameter space, the allowed mass matrix has no clear structure: All the elements are 
of the same order and could be taken as random numbers. This fact, and the existence of 
two large mixings, motivates the anarchy approach 48| where all parameters of the neutrino 



mass matrix are allowed to take random values. One then calculates the probability that 
012, ^23 and satisfy the experimental results. There can be variations to this approach 
depending on whether the anarchy is assumed to be in the high-scale sector of the theory 
(such as in the right handed neutrino mass matrix in the seesaw models described below) 
or in low-scale sector. One feature of these models is that in general they predict "large" 
^13 (closer to 0.1 or higher) for a large domain of parameters. One can also consider partial 
"anarchy" - randomness of some parameters on the top of a dominant structure determined 
by a certain symmetry. Questions for this approach include: which kind of physics leads to 
the anarchy and why does anarchy manifest itself for lepton mixings and neutrino masses 
only? To some extend, studying anarchy can be considered as a test of complexity behind 
neutrino masses and mixings. In fact, several comparable contributions to the neutrino mass 
matrix may exist, each having a simple structure and obeying a certain symmetry, but the 
totality of it having the anarchy effect. 

F. RenormaHzation group equation effects 

If the underlying theory is formulated at some high-energy scale M, e.g., much above the 
electroweak scale, one needs to use the renormalization group equations (RGE) to extrap- 
olate the mass matrix from the low-energy scale, where it has been reconstructed, to the 
high-scale M. In general, renormalization can change the structure of the Yukawa coupling 
matrix. Thus, to uncover the mechanism of mass generation one needs to calculate RGE 



corrections 49|]. In fact, some features of the mass matrix and observables at low scale can 



26 



be due to RGE effects. 

Suppose neutrinos are Majorana particles and their masses are generated at the EW scale 
by the operator eq. ([T]) (For renormalization of masses of the Dirac neutrinos see Ref. 50|.) 
Then, in general, two RGE effects should be taken into account: 

1) . renormalization of the operator in eq.([T]) from the low scale up to the scale at which 
it is formed, and 

2) . renormalization between and above the scale of the operator formation 5ll- In 



general, different terms of the operator in eq.([T]) are formed at different scales separated 
by many orders of magnitude rather than at the unique scale. This happens, e.g., in the 
case of the type I seesaw mechanism with strongly hierarchical masses of the right handed 
neutrinos: Mi -C M2 -C M3. The underlying physics is formulated at M > M3. In this case 
one should take into account threshold effects - a different renormalization group running 
between the masses and also above the seesaw scales (see sect. 6). 

In this section we consider the RGE effects below the scale of operator eq.([l]) formation. 
In the seesaw version with a hierarchy of the RH neutrino masses, that would correspond 
to /i < Ml, where /i is the running scale. 

In the SM as well as Minimal Supersymmetric Standard Model (MSSM), treatment of the 
RGE effects on the mass matrix in the flavor basis is simpler than RGE effects on observables 
- angles or masses. The observables can be found after renormalization by diagonalizing the 
obtained mass matrix (matrix of the Yukawa couplings). 



The RG equation for the effective mass matrix has a very transparent structure 49|, l52 1 

^ = CiY.WiM, + M,CiY,Wi + aM., (46) 

where t = (l/lGTr^) log(/i//io), Q = -3/2 in the SM and Q = 1 in MSSM. The first two 
flavor dependent terms correspond to the neutrino wave function renormalization due to 
Yukawa couplings of the charged leptons, the last term is the flavor independent renormal- 
ization due to gauge couphngs and also Yukawa coupling renormalization of the Higgs field 
wave function. 
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1. Renormalization of the neutrino mass matrix. 



In lowest order the gauge couplings produce the overall renormalization of the mass 
matrix only and do not change its flavor structure. (This is not true for threshold corrections 



because couplings of different right-handed neutrinos are flavor dependent [5l|.) In contrast, 
the Yukawa interactions modify the flavor structure of the mass matrices. In the flavor basis, 
Yukawa corrections do not generate the off-diagonal elements of the charged lepton mass 
matrix in the SM and MSSM. This matrix remains diagonal and therefore the correction 
does not change the flavor basis. On the contrary, RGE corrections change structure of the 
non-diagonal neutrino mass matrix. 

To understand the RGE effects we consider the one-loop corrections and neglect all the 
Yukawa couplings except the tau lepton one, Yr = rrir/vd, where Vd is the VEV that generates 
masses of down fermions. Essentially, the RGE effects are reduced to the wave-function 
renormalization and can be written as 

MM = IMR{fi)Mu{mz)R{fi). (47) 

Here Ic is flavor independent renormalization factor and 

R ^ diag{l, 1, _ 1 = log (48) 

The size of the effect is different in the SM and MSSM: In the SM Vd = v = 265 GeV and 
Yr = 9.5 X 10^^ at the EW scale, so the corrections are very small: For fi = 10^^ GeV 
we obtain [Z^ — 1)sm — 10^^. The effect is strongly enhanced in MSSM with large tan/3, 
where Vd ~ v/ tan/3, so that {2^ — 1)mssm ~ (^t — l)5Aftan^/3. For tan/? = 50, we obtain 
Zr-lr^ 0.03. 

Corrections appear as factors multiplying the bare values of the matrix elements. Thus, 
the zero elements will remain zeros [sSj. This allows us to draw important conclusion: The 
RGE effects (at least in the SM and MSSM) do not change the structure of the mass matrix 
significantly. Therefore the mass matrix reconstructed at low energies will have nearly the 
same form at high scales (before threshold corrections are included). This is not true if some 



new interactions exist above the EW scale (see e.g. 54|). 

In contrast, effect of corrections on the observables - angle and mass differences - can be 
strong for particular forms of the zero order mass matrix. 
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2. Renormalization of observables. 



The strongest effect on the observables is in the case of the quasi-degenerate mass spec- 



trum 
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i, 156|, 157J. Indeed, the corrections are proportional to the absolute mass scale mo: 
6m = tuq^Zj. — 1). It generates the mass squared difference Am^ ^ 2mo5m = 2mQ(Z^ — 1). 
Effect increases as square of the overall mass scale. For mo = 0.3 eV and {Zr — 1) ~ 10^^ 



it can give the solar mass split 



571]. Furthermore, the mixing angles depend on mass dif- 



ferences whereas corrections are proportional to the absolute values of the elements so that 



the relative corrections to the mixing get enhanced 



56| by A^^ oc m/Am oc mg/Am^. Es- 



sentially, the enhancement of mixing occurs when neutrinos become even more degenerate 



at low energies. In the case of normal 
Let us summarize possible effects 



hierarchy, however, the effect of the RGE's is small. 
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56 



57 



58 



59, 



60|. 



1). The angle 612 can undergo the strongest renormalization since it is associated to 
the smallest mass split. Some important dependences and results can be traced from the 
approximate analytical expression for running 58| : 

|mi + 17126'^^'^^'^ 



612 ~ -A sin 26112 sin^ ^23- 



Am|i 



+ 0(6 



13 J 



(49) 



where A = CiY^ /32tt'^ and 0i2 = 02 — 0i- Notice that other parameters in this formula, and 
especially Ami2, also run and their dependence on renormalization scale /i should be taken 
into account. As a result the dependence of 612 on the logfi turns out to be nonlinear. 

According to eq. (H9!) for 6*13 = 0, with increase of /i the angle 6^12 decreases in the 
MSSM and increases in SM. For the degenerate spectrum the enhancement factor can reach 
Ami/ A 

m'^2- Maximal enhancement corresponds to zero relative phase, 0i2 — 0, and running 
is suppressed for the opposite phases. 

For the degenerate spectrum the angle 612 can run practically to zero. This means that 
the large 1-2 mixing at low energies may have the radiative origin being small at, e.g., the 



GUT scale 
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56|. Another interesting possibility is that at Mqut equality of the quark and 
leptonic 1-2 mixings, 612 = 6'c7, is realized. So, the difference of quark and lepton mixings is 
related via the RGE running to degeneracy of the neutrino spectrum. 

In MSSM the angle 612 can increase with due to effect of non-zero ^13 provided that 
012 = vr, which equality ensures that the effect of the main term in eq. (jUj) is suppressed. 
2). Evolution of 1-3 mixing associated to the larger mass split is weaker and nearly linear 
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in log/i. It can be approximated as 

7TZ3 

6'i3 ^ A sin 26^12 sin 2623 n — 5- ["^1 cos(0i — 6) — 1712 cos(02 — 5) 

Amis 

-rAmscos^] + C(^i3). (50) 

The enhancement factor for the degenerate spectrum is mg/ Am^g and the strongest evolution 
is when phases 0i and 02 are different. For equal phases running is suppressed by an 
additional factor r^. In the case of normal mass hierarchy ^13 decreases with log/x. 

The main term in eg. (1501) does not depend on 613 at all, and therefore it evolves to nonzero 
value even if 613 = at some high energy scale. Consequently at low scales non-zero ^13 may 
have purely radiative origin. For instance, one can get ^13 ~ 8°, i.e. at the level of present 
experimental bound, if it is zero at Mqut 



If the mass hierarchy is inverted ^13 increases with /i and at the GUT scale it could be 
larger than 9c- The RGE effect is strongly suppresses when is small. 
3). Running of the 2-3 mixing can be described approximately by 



\ 

023 ^ -Asm2e23-^^[sm^ Oulmie'^' + rus]'^ 
Amis' 

+ cos^ ^121^26^2 + map] + 0(^13). (51) 

As in the previous case the enhancement factor for the degenerate spectrum is mg/Am23 
and running is suppressed if 0i = 02 = tt. In MSSM with increase of log /i the angle 623 



decreases and can be as small as (20 — 30 



radiative enhancement of the mixing 49 , 
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at the GUT scale. So, one can obtain the 



56|: 623 is small (similar to 9c) at high energies 



and it reaches ~ 45° at low energies. Here again the large lepton mixing is related to the 
neutrino mass degeneracy. 

The RGE should lead to deviation of the 2-3 mixing from maximal when running to small 
scales if it is maximal at high scale: D23 = \{Zr — 1) [goI], though for normal hierarchy the 
deviation is below 1°. 

Finally let us consider the renormalization of the 1-2 mass split: 

Ami2 ^ aAml^ - AA[2 sin^ ^23(^2 cos^ 612 - ml sin^ 612) + Oi^dvi)], (52) 

where the first term is the overall renormalization of all masses due to the gauge radiative 
corrections and also renormalization of the Higgs boson wave function. The second term 
can dominate for the degenerate spectrum. Depending on parameters the renormalization 
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can enhance splitting up to the atmospheric one or suppress it down to zero. So, zero 
spli t at some high scales and radiative origin of the 1-2 split at low scales can be realized 



571]. Another possibility is that at high scales all the mass splits are of the same order and 



the hierarchy of splits at low scales is produced by radiative corrections in the case of the 
degenerate spectrum. 

As is clear from our discussion the role of the RGE effects depends on a number of 
unknowns: possible extensions of the SM like two higgs doublet model, MSSM, (ii) on 
value of tan/3, (iii) on type of spectrum (degenerate, hierarchical), (iv) on CP violating 
phases. Depending on these unknowns the effects can vary from negligible to dominant, 
thus explaining main features of the neutrino mass spectrum and mixing. Apparently many 
uncertainties related to RGE effects will disappear if it turns out that the neutrino mass 
spectrum is hierarchical. Even in this case the corrections can be larger than accuracy of 
future measurements of the neutrino parameters. Strong effects are expected also from the 



"threshold" corrections 



5l|. 



G. Searching for the symmetry basis 

The structure of the mass matrix and its symmetries depend on the basis. In general, 
symmetry basis can differ substantially from the flavor basis considered so far. Therefore 
identification of the symmetry basis is crucial for uncovering the underlying physics. Un- 
fortunately, there is no clear guideline on how to look for this basis. One can perform a 
continuous change of the basis, searching for situations where both neutrino and charged 
lepton mass matrices have certain common symmetries. Some hints can be obtained from 
explicit models. 

It is not impossible for the symmetry basis to coincide with the flavor basis, and models 
where they do coincide exist. One can expect that in this case symmetries associated to 
some combinations of Lg, L^, play an important role. 

The symmetry basis may not coincide with but may be close to, the flavor basis. They can 
differ by rotation through an angle of the order of the Cabibbo angle, 6c ~ \J^^j./^t- The 
strong hierarchy of masses of charged leptons favors this possibility. However, the symmetry 
basis can strongly differ from the flavor basis. In some models, as a consequence of symmetry, 
the neutrino mass matrix is diagonal and maximal mixing comes from diagonalization of the 
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charged lepton mass matrix. Further studies in this direction are necessary. 



IV. NEUTRINOS AND NEW SYMMETRIES OF NATURE 

The most striking and unexpected outcome of the bottom-up approach is an indication 
of particular symmetries in the neutrino sector, namely, symmetries of the lepton mixing 
matrix and the neutrino mass matrix in the flavor basis. It is strange that symmetry is 
associated somehow with neutrinos and does not show up in other sectors of theory. Several 
observations testify to such a "neutrino" symmetries: 

- Maximal or close to maximal 2-3 mixing; 

- Zero or very small 1-3 mixing; 

- Special values of 1-2 mixing; 

- Hierarchy of mass-squared differences; 

- Quasi-degenerate mass spectrum, if established. 

As far as the last item is concerned, independent confirmation of the Heidelberg-Moscow 
positive result is needed. Note that in physics, large or maximal mixing is often related 
to degeneracy. THus , possibility of the degenerate spectrum does not seem implausible. 
Some of the features listed above can originate from the same underlying symmetry. 



A. fu — i^T symmetry 



Both maximal 2-3 mixing and zero 1-3 mixing indicate the same underlying symmetry 
and therefore deserve special attention. They are consequences of the — permutation 
symmetry of the neutrino mass matrix in the flavor basis 40|] . The general form of such a 
matrix is 

( A B b\ 
BCD 
B D C 



M 



V 



(53) 



This symmetry can be a part of discrete 5*3, or D4 groups and also can be embedded into 
certain continuous symmetries. The matrices in eq. fl36l) for c = 1 and a = b, and in eq. (HO|) 
for X = 2/ and c = s are special cases of the general matrix in eq. (l53|) . 
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At first, this might seen problematic because this symmetry cannot be extended to the 
charged lepton sector. To see this, note that in the fiavor basis for the charged leptons 
we have zero off diagonal elements, Di = Bi = 0, and therefore the symmetry implies 
Ml = diag{Ai,Ci,Ci) which contradicts <^ m^. If Di 7^ 0, implying that the symmetry 
basis does not coincide with the flavor basis, one can get the required mass hierarchy. 
However, in this case the charged lepton mass matrix also produces maximal mixing rotation. 
Neutrino and charged lepton rotations cancel leading to zero lepton mixing. 

Apparently the z/^ — Uj. symmetry should be broken. One way to resolve the problem is to 
couple the Higgs bosons, which violate yU— r symmetry spontaneously, weakly with neutrinos, 
but strongly - with the charged leptons. Examples where this is achieved using simple 



auxiliary symmetries such as Z2, have been discussed in the literature [61|- A difference 
between charged leptons and neutrinos appears because the right handed components Ir 
and uji have different transformation properties under Z2. This is not unnatural in the 
context of supersymmetric (SUSY) theories in which charged leptons get mass from the 
down Higgs, whereas the neutrinos get mass from the up Higgs doublet. Such models can 
be embedded into SUSY SU(5) GUTs \q^. 

Other possibilities are to introduce the symmetry basis which differs from the flavor 
basis (in this case the symmetry will be the 2-3 permutation symmetry), or to use other 
(approximate) symmetries which in the flavor basis are reduced to — permutation. 
Small breakings of these symmetries would manifest themselves in the appearance of a 
small but nonzero ^13. The question of course is what is small? It is reasonable to assume 
that ^^13 ~ ta (or smaller values of 1-3 mixing) indicates an underlying z/^ — symmetry. 
However, if Q\2 ~ a/^, no conclusion about this symmetry can be drawn because there are 
many examples in which larger values of 6*13 are possible. 

Can this symmetry be extended to the quark sector? Ref. [g^ argues that in fact 



smallness of the 2-3 quark mixing, V^h <^ yms/mfe, can also be a consequence of this 
symmetry. 

There are two shortcomings of the discussed symmetry: 

1) . It does not determine masses: Symmetry flxes general form of the mass matrix 
(equalities of certain matrix elements) and not masses which are given by the values of the 
elements. 

2) . Symmetry does not determine the 1-2 sector. Thus, one needs to use some more 
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extended symmetries which involve all three generations. One such symmetry that is widely 
studied is the ^4 symmetry |45|. Other possibilities are: 5*3 |6J], Z4 [65|, Da [66 |. 



B. A symmetry example, A4 



An interesting class o: 



of four elements 
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models is based on the A4 symmetry group of even permutations 



68| . It is the symmetry of the tetrahedron and has the irreducible 
representations 3, 1, 1' and 1". The products of representations 3x3 = 3 + 3 + 1 + 1' + !" 
and also 1' x 1" ~ 1 both contain the invariant 1. This allows one to introduce the Yukawa 
couplings with special flavor structure. Furthermore, it is the existence of three different 
singlet representations that leads to substantial freedom to reproduce the observed pattern 
of masses and mixings. Notice also that A4 is subgroup of 50(3). 

In all the models proposed so far, three lepton doublets form the triplet of A4: Li = 
(z/j, li) ~ 3, where i = 1,2, 3. The right handed components of the charged leptons, 1^, neu- 
trinos and Higgs doublets transform in different ways as either 3, or as 1, 1', 1", depending 
on the model. Essentially, the large (maximal) mixing originates from the fact that //j and 
have different A4 transformation properties. A disadvantage of the model is that it also 
requires the introduction of new Higgs multiplets, and often new heavy leptons as well as 
quarks which are generic features of most symmetry approaches. 

In A4 models one should introduce Higgs fields with non-trivial A^ transformation prop- 
erties, that is in representations 3 and 1, 1', 1". One can ascribe these properties to the SU2 
Higgs doublets, in which case 6 such doublets are required. Alternatively, one can keep SM 
Higgs doublet to be singlet of A4 but introduce new SU(2) singlet scalar fields which form 
non-trivial representations of ^4 (in the spirit of Froggatt-Nielsen approach). The latter 



however requires introduction of non-renormalizable operators (see e.g. 69|) or explicitly 
new heavy leptons and quarks [67]. 

There are different versions of the ^4 models. By appropriate choice of the Higgs fields 
and their VEVs and/or right handed neutrino couplings, one can obtain the tri-bimaximal 
mixing. The models constructed are based on the fact that tri-bimaximal mixing is given 
by the product of the trimaximal ( "magic" ) rotation Utm and maximal 1-3 rotation: 

Utbm = UtmU"^^. (54) 
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FIG. 6: Generation of the lepton masses and mixing in two models based on A4: a), model 



and b). model 



a 



7l[ |. Lepton and Higgs multiplets are in shadowed boxes and circles. Numbers at 
the boxes and circles indicate the A4 representations of the corresponding particles. The lines show 
the Yukawa couplings or bare mass terms in the models. We indicate also places where mixing is 
generated. 
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(55) 



As an illustration, in fig. [6] we show schemes of generation of the neutrino and charge 
lepton masses which lead to tri-bimaximal mixing in two models based on 



In the model a) |70| which is certain modification of the early proposal 45|, l67| , the two 
Higgs doublets H and H' are invariant under A4, and the SM singlets x form A4 triplet. 
New heavy leptons E, have to be introduced. The model content and transformation 
properties of the fields have been arranged in such a way that in the symmetry basis the 
charged lepton mass matrix produces the Utm rotation, whereas the neutrino mass matrix 
produces maximal 1-3 mixing 68|. The latter requires also certain VEV alignment. 

In the model b) [TlJ] the SU2 Higgs doublets form the triplet and singlets of A^. 

This illustrates generic problems and complexity of realization of "neutrino symmetries" . 
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C. Real or accidental 



The main question in this approach is whether the symmetries hinted at by observations 
are simply accidental or backed by real physics, especially in view of the fact that the 
price one pays in terms of the number of Higgs fields and/or extra leptons for constructing a 
realistic model can be very heavy. Thus, a fruitful approach is to look for possible deviations 
from symmetries in data {e.g., small 613 along with correlation between 613 and 623 — 7r/4) 
and explore models that may give other tests, because if symmetries are not accidental, they 
have consequences of fundamental importance: new structures are predicted, unification 
path may differ substantially from what we are using now, etc.. The symmetries may hold 
the key to understanding the flavor problem. 

Another possibility is that the symmetries are not accidental but the underlying theory 
has not been found yet and observed symmetry relations are hints of a new sector in physics, 
e.g. flavor- universal mixing with new singlets may produce a symmetric contribution to the 
active neutrino mass matrix (sec. 8). 

The only way to establish that symmetry is not accidental is to flnd new consequences 
of the symmetry, i.e. to make predictions in the context of a certain model and to test the 
predictions in experiment. It is important to flnd not just one but several predictions (see 
discussion in sec. 5.E). 

Finally let us stress that the observational basis for the existence of symmetries (real or 
accidental) is not yet well established. As we described in sec. 2.D, still signiflcant deviation 
from maximal 2-3 mixing is possible and 1-3 mixing could be relatively large. So, further 
experimental measurements will be decisive. 

V. LEPTONS AND QUARKS 

Joint consideration of quarks and leptons and searches for possible relations between 
quark and lepton parameters are of fundamental importance, because this may (i) provide 
a unifled clue for understanding fermion masses and mixings and (ii) give more insight into 
the uniflcation of particles and forces in nature. 

Below we compare quark and lepton masses and mixings and discuss various ideas about 
possible relations of quarks and leptons such as (i) quark-lepton symmetry; (ii) quark-lepton 
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unification; (iii) quark-lepton universality; and (iv) quark- lepton complementarity (QLC). 



A. Comparing leptons and quarks 

There is a strong difference between the masses and mixings in the quark and lepton 
sectors (fig. [7]). The ratios of masses of neutrinos and the corresponding upper quarks are 
mijvrif. < 10~^°, m^/mt < 10"^^. Lepton mixings are large, quark mixings are small. The 
1-2 mixing is the largest for quarks, whereas 2-3 mixing is the largest lepton mixing. The 
only common feature is that the 1-3 mixing (mixing between the remote generations) is 
small in both cases. 

More careful consideration, however, reveals some interesting features: It seems the 1-2 
as well as 2-3 mixing angles in the quark and lepton sectors are complementary in the sense 



that they sum up to maximal mixing angle 
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0i2 + 9c^^, (56) 



^23 + ^6-^. (57) 



Although for various reasons it is difficult to expect exact equality in the above relations, 
one can say, qualitatively, that there is a certain correlation: 2-3 mixing in the lepton 
sector is close to maximal mixing because the corresponding quark mixing is small, the 1-2 
mixing deviates from maximal mixing substantially because 1-2 (Cabibbo) quark mixing is 
relatively large. For the 1-3 angles we do not see a simple connection, and apparently the 
quark relation ^13 ~ 612 x 623 does not work in the lepton sector. Below we explore the 



implications of the above equalities known in the literature as QLC relations 721 ]. 

Comparing the ratio of neutrino masses in eq.( l23i) with ratios for charged leptons and 
quarks (at niz scale), 

m.j m, rrir , . 

— ^i = 0.06, — = 0.02- 0.03, — = 0.005 58) 

mi 

one concludes that the neutrino hierarchy (if it exists at all) is the weakest one. This is 
consistent with possible mass-mixing correlation, so that large mixings are associated with 



a weak hierarchy: Jmi/rrij ~ 6ij 
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FIG. 7: Mass hierarchies of quarks and leptons. The mass of the heaviest fermion of a given type 
is taken to be one. 



It is also intriguing that 



/ Hhi ~ gin 6*^2 ~ sin 6. 



c- 



(59) 



Does this perhaps indicate that the Cabibbo angle is the universal flavor parameter for 
both quark and lepton physics? In a class of grand unified models to be discussed below, 
the Cabibbo angle, indeed, becomes the key parameter of the neutrino mass matrix and 
describes the neutrino masses as well as mixings. 

In fig. [7] we show the mass ratios for three generations. The strongest hierarchy and 
geometric relation m„ x ~ exist for the upper quarks. Apart from that no simple 
relations show up. Furthermore, it looks like the observed pattern is an interplay of some 
regularities - flavor alignment and randomness - "anarchy" . Below we explore the possible 
meaning behind this picture. 



B. Quark-Iepton symmetry 

There are good reasons to suspect that quarks and leptons may be two different mani- 
festations of the same form of matter. The first hint arises from the observed similarities 
between weak-interaction properties of quarks and leptons. Each quark has its own coun- 
terpart in the leptonic sector, which has the same weak isospin properties: Ul corresponds 
to z//,, (11 - ei, etc.. This is generally known as quark-lepton symmetry, and although it 
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manifests itself only in the left-handed helicity sector of quarks and leptons, it is often con- 
sidered as a hint of further unification among these two very different kinds of matter. It 
can be extended to the right-handed sector, when leptons are treated as the fourth color {tS] 
following Pati-Salam SU{A)c unification symmetry. 

The second hint comes from the attractive hypothesis of grand unification of matter and 
forces which argues that at very short distances, all forces and all matter unify. In GUTs, 
quarks and leptons form multiplets of the extended gauge group. The most appealing such 
group is SO(IO) where all known components of quarks and leptons including the right- 



handed neutrinos fit into the unique 16-plet spinor multiplet [761]. It is difficult to believe 
that all these features are accidental. 

However, at first sight different patterns of masses and mixing strongly break the quark- 
lepton symmetry. Furthermore, if particular symmetries are found to exist only in the 
leptonic sector, they may indicate that quarks and leptons are fundamentally different. 

C. Quark- lepton universality 

In spite of the strong difference between the masses and mixings of quarks and leptons 
we still can speak about the approximate quark-lepton symmetry or even universality. The 
universality is realized in terms of mass matrices or matrices of the Yukawa couplings and 
not in terms of observables - mass ratios and mixing angles. The point is that very similar 
mass matrices can lead to substantially different mixing angles and masses (eigenvalues) if 
the matrices are nearly singular (approximately equal matrices of rank-1) jr^]. The singular 
matrices are "unstable" in the sense that small perturbations can lead to strong variations 
of mass ratios and (in the context of the seesaw) mixing angles. The well-known examples 



of singular matrices are the "democratic" mass matrix 78| and the matrix with only one 
non-zero element m^^. 

Let us assume that in the zeroth order of some approximation all fermion (quark and 
lepton) mass matrices are equal to the same universal singular matrix Yq: 

Y^ = Y, = Yd = Ym = Yi = Fo- (60) 

Here Yd is the Dirac-type neutrino Yukawa matrix. The Majorana type matrix for the RH 
neutrinos, Ym-, can in general differ from the others being, e.g., as Ym ~ Y^. 
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As an important example we can take 

/ A3 A2 \ 

Yo-- 



A3 A2 A 

l2 



, A~sin^c~0.2-0.3. (61) 

VA^ A 1 J 

This matrix has only one non-zero eigenvalue and no physical mixing appears because ma- 
trices of all fermions are diagonalized by the same rotation. In this respect all singular 
matrices are equivalent (up to basis definition) until corrections are introduced - corrections 
break the equivalence. 

Let us introduce perturbations of matrix structure, e, in the following form 

y/ = r,°(l + ef^.), f = u,d,e,u,N, (62) 
where Y^^^ is the element of the original singular matrix. This form can be justified, e.g., in 



the context of the Froggatt-Nielsen mechanism 



79| . It turns out that small perturbations 



e < 0.25 are enough to reproduce all large observed differences in mass hierarchies and 
mixings of quarks and leptons [t^]- 

The smallness of neutrino masses is explained by the seesaw mechanism. A nearly singu- 
lar matrix of the right-handed neutrinos, which appears in the denominator of the seesaw 
formula leads to enhancement of lepton mixing and can fiip the sign of the mixing angle 
which comes from diagonalization of the neutrino mass matrix. Thus,the angles from the 
charged leptons, Ue, and neutrinos, U^, sum up, whereas in quark sector, mixing angles from 
up- and down-quark mass matrices subtract leading to small quark mixing. 

Note that the different mass hierarchies of the up and down quarks (as well as charged 
leptons) may testify that two different universal matrices should be introduced for fermions 
with the third projection of weak isospin 1/2 and -1/2. This can also be related to existence 
of two different Higgs doublets giving masses to those components (as in SUSY theories). 

Instead of mass matrices one can consider the universality of rotation (mixing) matrices 
Vu which diagonalize mass matrices of all fermions in certain basis. The model can be 
arranged in such a way tha t in the lowest order Uckm = Vu^u = whereas Upmns = VijVu 



contains large mixings 80 1 



In connection to the quark-lepton universality one can consider the following working 
hypothesis: 
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1) . No particular "neutrino" symmetry exists, and in general one expects some deviation 
of the 2-3 mixing from maximal as well as non-zero 1-3 mixing. Nearly maximal 2 -3 mixing 
would be accidental in this case. Instead, some family symmetry is realized, which ensures 
universality of mass matrices and their particular structure. 

2) . The seesaw mechanism with the scale of riht-handed neutrino masses M ~ (10''— 10^^) 
GeV explains the smallness of neutrino mass. 

3) . The quark-lepton unification or grand unification are realized in some form, e.g., 
^O(IO). 

4) . The quark-lepton symmetry is weakly broken with still some observable consequences 
such as TTih ^ ttIt-. 

5) . Large lepton mixing is a consequence of the seesaw mechanism - seesaw enhancement 
of lepton mixing (special structure of the right-handed neutrino mass matrix), and/or due 
to contribution from the type II seesaw (which we consider in sect. 7.1). Flavor symmetry 
and/or physics of extra dimensions could determine this special structure. 



D. Quark-lepton complementarity 



As noted in eqs. fl56|) and fl57l) . the latest determination of solar mixing angle gives 

e^2 + 0c = 46.7° ± 2.4° (la) (63) 



which is consistent with maximal mixing angle within la (fig. [T]). The fact that the ap- 
proximate complementarity is als o fu lfilled for the 2-3 mixings hints at more serious reasons 
than just numerical coincidence 190|]. A possibility that the lepton mixing responcible for 



solar neutrino conversion equals maximal mixing minus 6c was first proposed in 82|, and 
corrections of the bimaximal mixing by the CKM type rotations discussed in 



If not accidental, the quark-lepton complementarity would require certain modification 



84l |. It implies the existence of some 



of the picture described in the previous section 
additional structure in the leptonic (or quark?) sector, which generates bi-maximal mixing. 
In this sense it might indicate a fundamental difference between leptons and quarks. How- 
ever, there should also be the quark- lepton unification or symmetry, which communicates 
the quark mixing to the lepton sector. A general scheme could be that 



"lepton mixing = bimaximal mixing — CKM". 



(64) 
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(Another option is "CKM = bimaximal - PMNS", which may have different imphcations) . 

There are a number of non-trivial conditions that must be met for the exact QLC relation 
to be realized [t^ . 

1) . The matrices f/]^ and Ui^^^'^ should be multiplied in the following order: 

UpMNs ^ U}U, = ...U^,...Ur2Ug''''^ (65) 

(last two matrices can be permuted). Different orders of rotations lead to corrections to the 
exact QLC relation. 

2) . A matrix with CP violating phases should not appear between [/{^'^^^ and U]^ or 



the CP violating phase in this matrix should be small |74j. |85|. 

3). The RGE effects should be small because presumably the quark-lepton symmetry 
which leads to the QLC relation is realized at high mass scales. 

Let us first describe two possible (to some extent extreme) scenarios of eq. ( IMl) that differ 
by the origin of the bi-maximal mixing and lead to different predictions. 

(1). QLC-1: In the symmetry basis maximal mixing is generated by the neutrino mass 
matrix: f/^ = Ubm] it can be produced by the seesaw mechanism. The charged lepton 
mass matrix gives the CKM rotation Ui = Ul;j^f^j, as a consequence of the quark-lepton 
symmetry: nii = rrid- In this case the order of matrices in eq. fl65l) is not realized; U^^^ 
should be permuted with and consequently the QLC relation is modified: sin0i2 = 
(1/V2) cos - 0.5 sin 

or 

sin ^12 ^sin(7r/4-0c) + O.5sin0c(V2-l)- (66) 

Numerically we find sin^6'i2 = 0.331 which is practically indistinguishable from the tri- 
bimaximal mixing value. The predicted 1-3 mixing, sin ^13 = sinOc / is close to the 
upper experimental bound (fig. [3]). Combining this with expression for 1-2 mixing we get 



an interesting relation 612 ~ 7r/4 — ^^13 (85|. 



2). QLC-2: In the symmetry basis maximal 1-2 mixing originates from the charged lepton 
mass matrix, Ui = Uhm, and the CKM, Ui, = U^j^j^^, appears from the neutrino mass matrix 
owing to the quark-lepton symmetry: rriD rriu (assuming also that in the context of the 
see-saw the right-handed neutrino mass matrix does not influence mixing, e.g., owing to 
factorization). In this case the QLC relation is satisfied precisely: sin ^12 = sin(7r/4 — 6c). 
The 1-3 mixing is very small - of the order sin^ \Vcb\ . 
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According to fig. [T] the best fit experimental value of 612 is in between the QLC-1 and 
QLC-2 predictions and further measurements of the angle with accuracy A612 ~ 1° are 
required to disentangle the scenarios. 

Other possibilities exist as well. For instance, one maximal mixing may come from the 
neutrino mass matrix in the symmetry basis and another from the charged lepton mass 
matrix. 

There are two main issues related to the QLC relation: 

- origin of the bi-maximal mixing; 

- mechanism of propagation of the CKM mixing from the quark to lepton sector. 

The main challenge here is that the required quark-lepton symmetry is broken. In par- 
ticular, the leptonic mass ratio rrie/m^ = 0.0047 is much smaller than the quark ratio 
uid/ms = 0.04 — 0.06; also the masses of the muon and s-quark are very different at the 
GUT scale. Precise QLC relation may imply that 

- the q-1 symmetry is actually weakly broken as we discussed in sec. 5.3; 

- the q-1 symmetry is very weakly broken for up quarks and neutrinos in a sense that for 
Dirac matrices ^ Md. Then CKM propagates via the up-sector; 

- the breaking affects mainly the masses and mass ratios but not mixings. 

Anyway, the mass matrices are different for quarks and leptons and propagation of 
the CKM mixing leads to corrections to the QLC relation at least of the order A612 ~ 
ecmd/rris ~ 0.5-1.0° [T^- 

Consider the case of QLC-1 (bimaximal mixing from neutrinos), where deviation of quark 
mixing from zero and lepton mixing from maximal follow from the down quarks and charged 
leptons. If the leptonic mass matrix has similar structure to the d-quark mass matrix with 
Gatto-Sartori-Tonin (GST) relation one would expect 61 ~ ^Juie/m^ ^ 6'c/3 and deviation 
from maximal mixing 0i/\/2 = l/3-\/2^c turns out to be too small [86]. There are several 
proposals to enhance the shift angle. In particular, the neutrino mass matrix can be modified 
a.s Aiu = Mlbm + SM, where ^Abm produces the bi-maximal mixing and 66 A4 leads to 



deviation [86|]. 6Ai can be due to the seesaw type-II contribution [87|]. However in this case 
connection to quark mixing is lost and the relation eq. fl56l) is simply accidental. Notice that 
the ratio of the mass squared differences, ta ~ sin6'c, so that the shift, 9i, can be related 
simply with generation of the solar mass split and therefore be of purely leptonic origin. 
In the context of quark-lepton symmetric models, the enhancement may have the group 
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theoretical origin. In 88|] for certain operators generating fermion masses the relation 6i = 
360/2 has been found, where factor 3/2 is the ratio of Clebsh-Gordon coefficients. 

The renormalization group effects on 1-2 mixing are in general small and furthermore 
they lead to increase of the angle 612 at low scales. The negative shift can be obtained from 
renormalization group effects in presence of the non-zero 1-3 mixing [74]. Also the threshold 
corrections due to some intermediate scale physics like low scale supersymmetry can produce 



the negative shift thus enhancing the deviation from maximal mixing 



891]. Strong shift can 



also be obtained from RGE effects between and above the seesaw scales related to the RH 



neutrino masses (see sec. 6.3) [51|. 



To avoid the additional 1/3 suppression of 61 one can abandon the GST-type relation for 
charged leptons. Then 61 ~ 6c would imply nearly singular character of the 1-2 leptonic 
submatrix. 

As remarked before, the quark-lepton symmetry can propagate 6c to lepton sector ex- 
actly if the neutrino mass matrix is the source of both bi-maximal mixing and the CKM 
rotations. The charged lepton and down quark mass matrices should be diagonal, and as a 
consequence of the quark-lepton symmetry, m„ = itld- The left rotations for these matrices 
give UcKM and the rest of the seesaw structure generates the bimaximal mixing. In this 
case, however, the GST-relation in the quark sector becomes accidental. If the bi-maximal 
mixing is generated by the charge leptons (lopsided scenario, see sec. 6.4.4) the QLC relation 
becomes precise [tsI- 

The role of CP-phases can be important in the q-1 relations [t^. CP violating phase 
in UcKM produces very small effect on QLC due to smallness of Vub- Also in this scenario 
the leptonic CP phase is very small. On the other hand appearance of the phase matrices 
in between U^^^^ and f/{^ will both modify the QLC and the leptonic Dirac phase, 5. 
Apparently, the relation between these two modifications should appear. In the QLC-1 
scenario insertion of the phase matrix Iphase = diag{e^°',e^^, 1), between two 1-2 rotations: 
U'-'^^^ IphaseUbm, Icads to the following change of the QLC relation [85|: 

6u^l-^cos{S-n). (67) 

So, the phase diminishes the shift, thus destroying the relation. Maximal shift required by 
QLC implies 6 ^ tc, that is, suppressed CP violation phase. 

There are few attempts to construct consistent quark-lepton model which reproduces the 
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QLC relation. The simplest possibility is the SU{2)l x SU{2)r x SU{£}c model which im- 
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The strategy 



plements the quark-lepton symmetry in the most straightforward way 
is to obtain (using an additional flavor symmetry) the neutrino mass matrix with inverted 
mass hierarchy which leads naturally to the bimaximal mixing (QLC-1 realization). The 
quark-lepton symmetry provides equality of the mass matrices Mi = M^, and consequently 
the same CKM type mixing in both sectors. The perturbations to the matrices, 6Mi and 
SMd, should be introduced which break the q-1 symmetry and correct the masses. (In js^ 
they are due to the non-renormalizable operators with new Higgs fields which transform as 
15 of SU{A)c)- These corrections however modify relation between 6i and 6c, and their 
equality is matter of tuning of continuous parameters. Another possibility [88] is to intro- 
duce the non-renormalizable operators which include couplings with Higgs in 4 of SU{4)c 
as well as singlet fiavon fields a la Froggatt-Nielsen. Selecting particular type of operators 
one can get inequality of matrices Mi and M^ already in the lowest order and enhance the 
leptonic angle: e.g. like 6i = 36c/ 2, as we have marked previously. The enhancement allows 
to reproduce the QLC relation eq. (|56ll almost precisely. 

Different approach to resolve the problem of decoupling of masses and mixing is to use 
non-abelian flavor symmetries [73]. Via minimization of the potential the symmetries lead 
to zero or to maximal (bi-maximal) mixing independently of the mass eigenvalues. 

The Cabibbo mixing can be transmitted to the lepton sector in a more complicated 
way (than via the quark-lepton symmetry). In fact, sin6'c may turn out to be a generic 
parameter of the theory of fermion masses - the "quantum" of flavor physics, and therefore 
to appear in various places: mass ratios, mixing angles. The relation in eq. (!59l) favors of 
this possibility. However, the same relation in eq. (l59|) may suggest that the QLC relation is 
accidental. Indeed, it can be written and interpreted as pure leptonic relation 



Oi2 + 6^r = J, tan V = J—- (68) 

This relation may be even more difficult to realize in models. 

Following an idea that A ~ sin 6c is the "quantum" of the flavor physics one can consider 
the Cabibbo angle as an expansion parameter for mixing matrices. In zero approximation 
the quarks have unit mixing matrix: U^km = -^5 whereas leptons have Up^,j^g = Uhm 
or bi-large mixing matrix. The A size corrections can be included as Upmns = UlUtm 
or UpMNS = UhmU\- Interesting possibility (in a spirit of the QLC relation) is that 
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Ux = UckmW in the Wolfenstein parameterization [83|, l90|. In this case one gets uni- 
versal description (parameterization) of quark and lepton mixing matrices. This apparently 
reduces the number of free parameters in the problem and also establishes various relations 
between mixing angles. 

In general on can take U\ as a matrix with all three A size rotations and study properties 
of the PMNS matrix obtained by insertion of the Ux in various places of the zero order 



structure 



that is, UlUbm, UbmUl or U^^UlU^^, etc. 



E. Empirical relations 

Establishing empirical relations between masses and mixings of fermions may give a clue 
to the underlying physics. The tri-bimaximal mixing scheme and QLC equality are examples 
of relations "between mixings without masses" . One should note, of course, the GST relation 



sin6'c ~ Vus ~ Jnid/fng (92|, and vn^jms = Jmu/rric which determine substantially the 



form of quark mass matrices, etc. 

A particularly intriguing such a relation is the Koide relation 
pole masses of charge leptons 
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95| between the 



^ mp + m.j + 2 , , 

Qi = , = - 69 



which is satisfied with accuracy 10 ^: 

^{pole) _ 9 /q +0.00002 (nr^\ 

— ^/o _o. 00001- V^i 



The Koide formula eq.(!69l) is interesting not only because of precision but also because it 
was obtained in the context of certain model and not as empirical relation. In fact, it was 
obtained in attempt to explain the relation (also rather precise) between the Cabibbo angle 
and lepton masses in the composite model of the quarks and leptons [94,]- It allowed to 
predict precise value of the tau-lepton mass. 

There are several properties of the relation eq. fl69|l which could have interesting implica- 



tions [961] 



(i) Varying masses one finds that the minimal value, Qmin = 1/3, corresponds to the 
degenerate spectrum and the maximal one, Qmax = 1, to the strongly hierarchical spectrum. 
So, the quantity Qi is a good measure of degeneracy of spectrum. The experimental value 
2/3 is exactly in between the two extremes. 
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(ii) The relation involves 3 generations explicitly. The mass of electron can not be ne- 
glected and therefore in the underlying theory rrie can not be considered as perturbation. In 
fact, the value 2/3 may be interpreted as 2/Nf, where Nf = 3 is the number of flavors. 

(iii) The formula is invariant under interchange of flavors e ^ /i, etc., and therefore implies 
5*3 (or wider) underlying symmetry. The value 2/3 may have certain group theoretical origin. 

(iv) Essentially eg. (1691) gives relation between the two mass hierarchies Ve = me/mr, 
= m^/rrir, and does not depend on the absolute scale of masses: 

0_ ^ + re + r, _2 

" + " 3- ^^'^ 

So, it can be realised for different sets of hierarchies. 



(v) The formula may have certain geometrical origin 97|, lOSj. Introducing vectors M = 
(1, 1, 1) and L = [y/m^, ^m^, y/rrir) we can rewrite it as 

1 . L-M , , 

3cos2 6'ml 1^1 1^1 

Apparently the experimental result corresponds to Qml = 45°. 

(vi) The relation has bilinear structure in ^/m which may imply that masses are bilinear 
of some other physical quantities: coupling constants or VEV's. In fact, the Koide relation 
is reproduced if 

rrii = mo{zi + ZqY, (73) 



where zq = yj^i^f/'^ and J^i^i = 0- Such a situation can be realized in the case of the 
radiative mechanism of mass generation: in one loop m oc F^, or in the seesaw mechanism 
m oc fifi'M~^. 

(vii) The quantit y Qi is not invariant under RGE running. At the - mass Qi{mz) ~ 



1.002QSp°''^ |99|,1iQ0|. Above mz the renormalization is negligible in the SM, and it can lead 
to further increase of Qi by about 0.7% in MSSM at Mr ~ 10^^ GeV and for tan /? = 50 |lOO |. 
So, the renormalization effect is much larger than the error bars in eq. flTOj) and therefore 
Qi deviates from 2/3 at high scales (already at the EW scale). This may indicate various 
things: the relation is accidental; the accuracy for the pole masses is accidental; physics 
responsible for the relation, and therefore the lepton masses, is at low scales. 

(viii) The relation eq. (!69|) is not universal: it can still be valid for the down quarks: 
Qd ~ 0.7 at mz, but it is certainly violated for the up quarks: Qu ~ 0.9. 
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(ix) Important aspect is that the mass relation does not depend on mixing. That is, 
physics of mass generation and that of mixing should decouple. Mixing can be included if the 
relation Q = 2/3 is considered, e.g., for the "pseudo-masses" introduced as = J2i Uiai^i, 
where Ul is the matrix of rota tion of the LH components which diagonalizes the mass matrix 
in the "symmetry" basis 101| . For charge leptons one should take = /. For quarks one 
can select the matrices so that universality = Q*^ = 2/3 is restored. 

Till now no realistic and consistent model for the Koide relation is constructed (see 

oop) 



96l | for review). Among interesting proposals one should mention the radiative (one 
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9J|; the seesaw mech anism 



95 



102|; 



mechanisms of charged leptons masses generation 
mechanism based on the democratic mass matrices and S3 symmetry |l03j |. An interesting 
possibility is that lepton masses are generated by bi-linear of VEV's of new scalar fields: 
rrii oc {(pi) {(pi) . Then as a consequence of symmetry of the scalar potential (6*3 and S'f/(3) 
symmetries have been considered), the VEV's have the property {(pi) Zi + zq in eq. (173|) 



96|. 



What about neutrinos? On the first sight, because of weaker mass hierarchy eq. (l22il . 
the neutrino masses do not satisfy the Koide relation. Depending on the unknown absolute 



mass scale one finds Qi^ = 0.33 — 0.60 
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lOOj], where the lower bound corresponds to 



the degenerate spectrum and the upper one to mi = 0. However, it was noticed recently. 



that the relation can be fulfilled provided that ^/m^ < and two others are positive 
So that for neutrinos the relation is 



104| 



Qu 



Till + 1712 + ^3 



2 
3' 



(74) 



Of course, the question is still why only one lepton z/i has negative root squared of mass. 
The relation (17^ together with the best fit values of the mass squared diffe renc es implies 



strongly hierarchical mass spectrum with prediction for the lightest mass 



(105|) 



104| (see also 



mi 



3.9-10-^ eV. 



(75) 



Two other masses m2 = 9.0 • 10^^ eV and = 5.1 ■ 10^^ eV are given essentially by the 
measured mass squared differences. 

Other proposals include the following: The universality can be restored if one uses the 



pseudo-masses 



10l[ |. Notice that since Uj^ = I for charged leptons, for neutrinos we have 
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Ul = UpMNS- Then from the condition Qu = 2/3 one finds for the allowed region of neutrino 
oscillation parameters: mi ~ (3 ± 1)10^^ eV, Ou > 35° and 623 > 50°. All neutrino masses 
are of the same order, and large lepton mixing is related to the absence of mass hierarchy a 
la the GST relation. 

Another proposal is to modify the Koide relation for the upper quarks and neutrinos 
without introduction of mixing [99|. Observing that Qu < 2/3 but Qu > 2/3 one can 
assume a kind of mass complementarity Qi + Qd = Qu + Qu- That would lead to the lightest 
neutrino mass mi ^ 10~^ eV. 

Notice that in these considerations smallness of neutrino mass and its possible Majorana 
nature have not been taken into account. Apparently, the presence of the Majorana mass 
matrix of the RH neutrinos in the context of seesaw mechanism can infiuence the implica- 
tions of the Koide relations for neutrinos. Alternatively one can imagine that mechanism 
responcible for smallness of the neutrino masses does not infiuence ratios of masses. 

The question: "real or accidental" is still open; and the lesson is that just one very precise 
prediction confirmed by very precise measurements may not be enough to verify theory. 



VI. SEESAW: THEORY AND APPLICATIONS 

As noted above, the see-saw mechanism is one of the simplest ways to understand the 
small neutrino masses. It has important implications and connections to a number of fun- 
damental issues which we discuss in this section. 

• What is the scale of Mji and what determines it? 

• Is there a natural reason for the existence of the right handed neutrinos, the main 
element of the see-saw? 

• Is the see-saw mechanism alone enough to explain all aspects of neutrino masses and 
mixings? 

• On a phenomenological level, what is the fiavor structure of the right handed neutrino 
sector? Can we determine it from, for example, purely low-energy neutrino observa- 
tions? 
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A. Right-handed neutrino masses and scale of seesaw 

The scale of the seesaw (type I) is related to the scale of right handed neutrino masses. 
Some idea about can be obtained from the naive estimation of masses for the third 
generation: 



where nit is the top quark mass, k{Mji) is the renormalization group factor of the D = 5 
operator. (Here we assume normal mass hierarchy.) It is this large scale that indicates that 
neutrino mass is related to new physics beyond that implied by the charged fermion masses. 
The scale in eq. (1761) is rather close to the GUT scale and in fact can be immediately related 
to the GUT scale. In this sense the smallness of the neutrino mass is the direct indication 



The situation is more complicated if one considers all three neutrinos and takes into 
account mixing among them. For instance, for the hierarchical mass spectrum of light 
neutrinos the heaviest RH neutrino mass is determines by the smallest neutrino mass mi 
and not m^: M3 ~ m^/nii. If mi < 10~^ eV, the mass M3 can be at the GUT scale or 
even bigger. There are a number of uncertainties and ambiguities in determining Mr. (i) 
we do not yet know the scale of light neutrino masses (which can change the estimation by 
approximately one order of magnitude); (ii) the Dirac masses of neutrinos are not known and 
must be assumed; (iii) mixing can strongly influence the masses of right-handed neutrinos; 
(iv) it is not clear yet that seesaw type I gives the main contribution to the neutrino mass. 
If it is subdominant, the masses of the right-handed neutrinos can be larger; (v) more than 
three singlet fermions (RH neutrinos) can be involved in the generation of the light neutrino 
masses. In this case (as we see below) the scale in eq. (1761) may turn out to be a phantom 
scale that does not correspond to any physical reality. 

The assumption that m^ oc m^ leads typically to a rather strong hierarchy of the RH 
neutrino masses: Mj oc m^j or even stronger. Their values can cover the interval (10^ — 10^^) 
GeV, although in some particular cases two masses can be quasi-degenerate. To get small 
masses for the usual active neutrinos, it is sufficient to have only two RH neutrinos, which 
means that the third one can be arbitrarily heavy: e.g., at the GUT or even Planck-mass 
scale. 




(76) 



of GUT. 
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B. Seesaw, B - L and Left-Right symmetries 



What is the physics content of this new scale? The seesaw scale can be identified as the 
scale of violation of certain symmetries. The fact that can be much smaller than the 
Planck scale is an indication in favor of this. It is therefore appropriate at this point to 
discuss possible origins of the RH neutrino masses. 

To answer this question, it is important to note the changes that occur in the standard 
model with the addition of one right handed neutrino per generation: The most obvious 
change is that it restores the quark-lepton symmetry. On a more fundamental level, it 
turns out that in the presence of three right-handed neutrinos, the symmetry B ~ L which 
was a global symmetry in the standard model becomes a gaugeable symmetry because the 
condition Tr{B—L)^ = is satisfied, implying that gauge anomalies cancel. The gauge group 
of weak i ntera ctions expands to become the left-right symmetric group SU{2)l x SU{2)ji x 



U{1)b-l |106| which is a subgroup of the SU{2)l x SU{2)fi x S'f/(4)c group introduced by 
Pati and Salam Tsl . This leads to a picture of the weak interactions that is fundamentally 
different from that envisaged in the standard model in that the weak interactions, like the 
strong and gravitational interactions, become parity conserving. Furthermore, in this theory, 
the electric charge formula becomes 107|: 

Q = hL + hR + ^-j^, (77) 

where each term has a physical meaning unlike the case of the standard model. When only 
the gauge symmetry SU{2)ji x U{1)b-l is broken down, one finds the relation AI^r = 
—A ("^2^)- This connects B — L breaking, i.e. A{B — L) ^ 0, to the breakdown of parity 
symmetry, AI^n 7^ 0. It also reveals the true meaning of the standard model hypercharge 
as 2" — ^3,R H — —■ 

To discuss the implications of these observations for the see-saw mechanism, note that in 
the first stage, the gauge symmetry is broken by the Higgs multiplets Al{3, 1, 2)©A/j;(l, 3, 2) 
to the standard model and in the second stage - by the bi-doublet 0(2, 2, 0). In the first stage, 
the right handed neutrino picks up a mass of order / < Ajj >= fvn, then (p produces the 
Dirac mass term. The presence of the coupling of the triplets with the bi-doublet 
leads to a shift of the minimum of the potential from Aj;^ = 0, so that this triplet acquires 
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FIG. 8: Feynman diagrams for a), type I and b). type II seesaw mechanisms. 

the so-called induced VEV of the size 

VL = (A° ) = ^ (78) 

from the diagram in figl8)D. As a consequence the mass matrix ([8]) is generated with the 
components 

rriL = fvL, mo = Yv^jk, Mr = Jvr. (79) 
The light neutrino mass matrix can then be written as 

M, = ^{Xf -Y^r^Y). (80) 



Note that vl and the see-saw type II term are suppressed by the same 



type I contribution, so that the overall see-saw suppression remains {g]. As a consequence 
of the left-right symmetry, the two contributions are partly correlated: Both depend on 
the same matrix /. An important problem is to reconstruct / from the low-energy data. 
In this connection, an interesting property of the formula in e a.f[80l) in the lowest order 
approximation (before RGE corrections) is the see-saw duality |l08l |: For any solution /, 
a dual solution / = M.v/vl — f exists. In the limit of very large right handed neutrino 
masses, the general seesaw formula for neutrino masses reduces to the triplet seesaw (type 
II) formula. 

Note that the scale of the right handed neutrino masses (or the seesaw scale) is now the 
scale of B-L breaking, which shows that it is not the Planck scale as would be the case if the 
left-right symmetry group were not considered. We see below that in the context of SO(IO) 
GUTs which embed the left-right model, the see-saw scale can indeed be the GUT scale, 
removing one arbitrariness in the description of neutrino masses. 
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'actor as the see-saw 



C. Seesaw and the RGE effects 

As we saw in sec. 6.1 the masses of RH neutrinos are in general substantially smaller than 
the GUT scale. Furthermore, typically they have an extremely large spread (often related to 
the quadratic mass hierarchy of quarks): from 10^ to 10^^ GeV. That determines important 
features of the RGE effects if the flavor physics (structure of the Yukawa couplings) is fixed 
at Mf > Mgut- In this case there are three different energy regions with different RGE 
behavior: 

1) . Region below the seesaw scales, /i < Mi where Mi is the mass of the lightest RH 
neutrino. The RGE effects in this region have been studied in sec. 3.6. 

2) . Region between the seesaw scales: Mi < fj, < M3, where M3 is the heaviest RH 
neutrino. 

3) . Region above the seesaw scales: /i > M3. 

In the regions 2) and 3) the key new feature is that some or all RH neutrinos are not 
decoupled and therefore the neutrino Yukawa couplings, Yi,, contribute to the running in 
addition to Fe- The couplings now run. The term C^YjY^, where C^, = 0.5 in SM and 
Cj, = 1 in MSSM, should be added to the RGE eq. (H^ . Also a should be modified. 

The couplings Yi, can be large - of the order 1 both in SM and MSSM independently of 
tan j3. Therefore in general the RGE effects due to Y^, are large. 

Another important feature is that the matrices Y^, and Y^ can not be made both diagonal, 
e.g., Yy is nondiagonal in the flavor basis where Ye is diagonal. This means that RGE 
running due to Fj, generates the flavor transitions and therefore leads to rotation of the 
flavor basis. This running can produce flavor mixing even if initial (at boundary) mixing 
matrix is proportional to unity [51]. 

In the region above the seesaw scales the running has similar features to those in the range 
(1). In particular, similar enhancement factors appear in the case of degenerate or partially 
degenerate spectra. Also CP-violating phases influence running substantially leading in 
certain cases to damping of the enhancement. The 1-2 angle can undergo the strongest 
renormalization. 

Let us note some interesting possibilities. For the degenerate spectrum and certain values 
of phases, running above the seesaw scales, (10^^ — 10^^) GeV can reduce 612 from 45° down 
to ~ 30°, thus explaining deviation of the lepton mixing from bi- maximal. It can correct the 
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QLC-1 relation reducing 612 at low energies. Running of other two angles is substantially 
weaker. Renormalization effects in two other region can be small, e.g., due to small tan/3. 

For the hierarchical mass spectrum (mi < 0.01 eV) the RGE induced changes of the 
mixing parameters are relatively small: A0i2 < (1 — 2)° (though it may be relevant for the 
QIC relation), Asin^^is < 3 ■ 10^^ and Asin^6'23 < 0.02. Mass squared difference Am^2 
can decrease due to running between and above seesaw scales by factor of 2 for partially 



degenerate spectrum, etc. 



5l|. 



In the region between seesaw scales one or two RH neutrinos decouple. The effective 
neutrino mass matrix has two different contributions: J^y = Airun + M.dec - the d=5 type 
term from the decoupled states, M.deci and the running term, M.run = {jj) M]^^ {jj)Yy{n) 
where where is the submatrix (3 x 2 or 3 x 1) of the Yukawa couplings for "undecoupled" 
RH neutrinos. In non-supersymmetric models these two contributions renormalize differently 
due to vertex corrections (in the SM case) to the d=5 operators. For the undecoupled RH 
neutrinos that corresponds to the converging box diagram with the RH neutrino propagator. 
The vertex corrections can change substantially the observables, e.g., leading to A^i2 ~ 10° 
even for the hierarchical spectrum 51|. 



D. Other realizations of seesaw 



If it turns out that the scale of B — L symmetry breakin g, Mr _/-,, is in the TeV range, as, 
for example, in a class of string models discussed rec ently 109| . the small neutrino masses 



can be understood by a double seesaw mechanism [110l | where, in addition to the right 
handed neutrino, A^, one postulates the existence of a singlet neutrino S. The symmetries 
of the model are assumed to be such that the Majorana mass of as well as the coupling 
of S to the lepton doublet are forbidden. We then have a neutrino mass matrix in the basis 
{u, N, S) of the form: 

^0 ^ 

M = rriD M . (81) 
^0 fi J 

For the case ^ <^ M Mb-l, this matrix has one light and two heavy quasi-degenerate 
states for each generation. The mass matrix of light neutrinos is given by 

My ~ mlM^-^i^M-^niD. (82) 
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There is a double suppression by the heavy mass compared with the usual seesaw mechanism, 
hence the name double seesaw. One important point here is that to keep fi ^ mo, one also 
needs some additional gauge symmetries, which often are a part of the string models. 

Another possibility which is motivated by the fact that the required masses of the RH 
neutrinos are at a somewhat smaller scale than the GUT scale, is that the RH neutrinos 
themselves get mass via a see-saw mechanism generated by and S. That would correspond 
to yu > M in eq.([8l]), so that 



1 T\jT 



Mr = -Mfi-'M 



(83) 



For fi ~ Mpi and M ~ Mguti that gives the required masses of the RH neutrinos. In 
particular, eq. fl83l) can produce a strong hierarchy of masses. The formula for the light 
masses is the same as in eq. (l82|) . 

If there is parity symmetry in models that implement the double seesaw mechanism, 
then the 13 and 31 entri es of the above neutrino mass matrix (IHTl) get filled by small see-saw 
suppressed masses |lll| . This leads to 

^ ml MV^ 
M = M , (84) 

^Me' /i J 

where e' ~ v^k/Vo and Vq is of the order the mass M in eq. (l84|) . The left handed neutrino 
mass in this case is given by 



= m^M'^-^fiM-^mD - {mo + mj^)e'. 



(85) 



The last contribution, linear in the Dirac masses, is termed t he see saw type III. There have 
been a few applications of this mechanism to model-building 112 |. 

There are also other variations on the seesaw theme for instance having two right handed 
neutrinos rather than three. Two RH neutrinos is the minimum number that will give a 
realistic spectrum for neutrinos after the see-saw mechanism. There are scheme s in which 
new symmetries beyond the standard model can realize such a possibility 47|, Ill3l |. For 
instance, if we supplement the standard model by a local SU {2)h horizontal symmetry that 
acts on the first two generations, then global anomaly freedom requires that there be only 
two right handed neutrinos transforming as a doublet under SU{2)h- Th is m odel leads to 



a 3 X 2 seesaw and has features similar to the two RH dominance models 



114 | 



55 



E. Seesaw and large lepton mixing 



1. See-Saw enhancement of mixing 



Can the same mechanism that explains the smallness of the neutrino mass, i.e., the see- 
saw also explain the large lepton mixing, so that eventually large mixing originates from 
zero neutrino charges and Majorana nature? The idea is that due to the (approximate) 
quark-lepton symmetry, or grand unification, the Dirac mass matrices of the quarks and 
leptons all have the same (or similar) structure: m^) ~ rriup ^ mi fn^own leading to zero 
(small) mixings in the first approximation. Owing to the non-diagonal mass matrix of RH 
neutrinos, Mfi, which has no analogue in the quark sector, the seesaw mechanism produces 
non-zero lepton mixing already in the lowest order. 

The problem with this scenario is the strong hierarchy of the quark and charged lepton 
masses. Indeed, taking the neutrino Dirac masses as m£, = diag{my_,mc,mt) in the spirit of 
grand unification, we find that for a generic Mr the see-saw type I formula produces strongly 
hierarchical mass matrix of light neutrinos with small mixings. The mixing becomes large 



only for special structure of M^j which 
Two different possibilities are 115 



compensates for the strong hierarchy in m^). 



116 | 



strong (nearly quadratic) hierarchy of the RH neutrino masses: Mm 
can be reproduced naturally by the double seesaw; and 

strongly off-diagonal (pseudo-Dirac) structure of Mji such as 



[m. 



tup J 



which 



M, 



R 



A 



(86) 



B 
\0 BO J 

which implies certain symmetry. Alternatively, the 12, 21 and 33 elements can be non- 
zero. An interesting consequence of these structures is that the pair of RH neutrinos 
turns out to be nearly degenerate, which can lead to the resonant leptogenesis. 

In the three-neutrino context both possibilities can be realized simultaneously, so that the 
pseudo-Dirac structure leads to maximal 2-3 mixing, whereas the strong hierarchy A B 
enhances the 1-2 mixing 117 1. 

There are several alternatives to the seesaw enhancement. 
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2. Large mixing from type II see-saw 



In general, the structure of the neutrino mass matrix generated by the type II (triplet) 
see-saw is not related to the structures of the matrices of other fermions and it can produce 
large mixing (see sect. 7). In some particular cases, however, the relations can appear 
leading to interesting consequences. 



3. Single RH neutrino dominance 

The large neutrino mixing and relatively strong mass hierarchy implied by the solar and 
atmospheric neutrino data can be reconciled if only one RH neutrino gives the dominant 



contribution to the see-saw |118[ |. (This leads to the 2-3 submatrix of my with nearly zero 



determinant.) There are two different realizations of this possibility. In one case the large 
mixing originates from the large mixing in the Dirac neutrino mass matrix mn: two left- 
handed neutrinos have nearly equal couplings to the dominating right-handed component. 
Suppose that (^,£1)23 ~ (?^_d)33 = "o^, (^d)i3 = Am (A ~ 0.2) and all other elements of rriD 
are much smaller. Then if only (M~^)33 is large in the inverted matrix, the see-saw will lead 
to a mass matrix with the dominant fi — t block. The mechanism can also be extended to 
enhance 1-2 in ixing. It requires the so-called sequential dominance related to the second RH 
neutrino 119|. 



In another version, the dominance is realized when two RH neutrinos are much heavier 
than the third (dominating) one and no large mixing in appears. This is equivalent 
to the strong mass hierarchy case of the see-saw enhancement mechanism. A realization 
requires {1710)22 ~ (^^0)23 ^ ("^0)33, and dominance of the (M~^)22 element. 

It may happen that the enhancement of the mixing is not related to the seesaw mechanism 
at all being, e.g., of the radiative origin. Let us consider the following possibility. 



4. Lopsided models 

Large lepton mixing in these models follows from the ch arged lepton mass matrix in the 



symmetry basis which should be left-right non-symmetric 120l |. This does not contradict 
the grand unification since in GUT models such as SU(5), the left-handed components of 
leptons are unified with the RH components of quarks: 5 = {d'^, d'^, d'^, I, v). Therefore, large 
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mixing of the left-handed leptonic components is accompanied by large mixing of the RH 
d-quarks which is unobservable. By introducing a Dirac mass matrix of the charged leptons 
in which the only large elements are {mi)^^ ~ (^^023 in the basis where neutrino mass matrix 
is nearly diagonal, one obtains the large 2-3 lepton mixing. This scenario can also be realized 
in 5*0 (10), if the symmetry is broken via SU{5). A double lopsided matrix for both large 
mixings (solar and atmospheric) is also possible. 



F. Screening of Dirac structure 

The quark - lepton symmetry manifests itself as a certain relation (similarity) between 
the Dirac mass matrices of quarks and leptons, and this is the origin of problems in explaning 
the strong difference of mixings and possible existence of neutrino symmetries. However, 
in the context of double seesaw me chanism t he D irac structure in the lepton sector can be 



completely eliminated - "screened" 1151 . ll2lL Il22| thus opening new possibilities. 



Indeed, the double seesaw mechanism leads to the light neutrino mass matrix given in 
eq. (!82|) . Suppose that due to a certain family symmetry or grand unification (which includes 
also new singlets S) the two Dirac mass matrices are proportional to each other: 

Md = K-'ruD, K = vew/Vgu- (87) 

In this case the Dirac matrices cancel in fl82|) and we obtain 

= K^iJ. (88) 

That is, the structure of the light neutrino mass matrix is determined directly by fi and 
does not depend on the Dirac mass matrix. Here the seesaw mechanism provides the scale 
of neutrino masses but not the flavor structure of the mass matrix. It can be shown that at 



least in the SUSY version, the radiative corrections do not destroy screening 121 |. 

The structure of the light neutrino mass matrix is given (up to small corrections) by n 
which can be related to some new physics at, e.g., the Planck scale. In particular, 

1) . fi can be the origin of neutrino symmetry; 

2) . /i oc / leads to a quasi-degenerate spectrum of light neutrinos; 

3) . /i can be the origin of bi-maximal or maximal mixing thus leading to the QLC relation 



8J] if the charged lepton mass matrix generates the CKM rotation (QLC-1). In general, 
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screening allows one to "automatically" reconcile the quark-lepton symmetry with the strong 
difference of mixings of leptons and quarks. 



G. Seesaw: tests and applications 

A major problem in neutrino physics is to find ways to test the proposed mechanisms 
and scenarios of neutrino mass generation. The seesaw scenarios are related to physics at 
very high energy scales which can not be achieved by the direct studies. Furthermore, it is 
practically impossib le to reconstruct the right handed neutrino mass matrix from the low 
energy ob serva bles 



neutrinos 



123[ | without additional assumptions like involvement of only two RH 

124j . etc.. 

The situation can change if the seesaw mechanism is embedded in into bigger picture so 
that one will be able to test the whole context. This will allow one to connect neutrinos with 
other phenomena and observables. Moreover, some parameters of the seesaw mechanism can 
be determined (see sect. 7). The hope is that in future on the basis of certain models we will 
be able to make predictions with very small uncertainties which can be tested in precision 
measurements. This will then provide a direct test of the model. 

In what follow we will describe briefly some connections of the seesaw with other phe- 
nomena which can help to check the mechanism. 



1. Testing see-saw in colliders 

See-saw mechanism combined with supersymmetry can provide another testable signa- 



ture in colliders 125l | by affecting the mass spectrum of sleptons at the weak scale. In 
simple MSSM scenarios, one assumes a universal scalar mass for all superpartners at the 
susy breaking scale (say Mqut oi Mp£ ). The weak scale masses are derived from this 
by renormalization group extrapolations, which are sensitive to the interactions at various 
scales. Above the see-saw scale, the neutrino Dirac couplings Y,y affect the e volu tion of both 



the left-handed slepton and sneutrino masses, as has been pointed out in |l25l |. They find 



effects which can be observable at the LHC Specifically, combination of slepton masses 

2A. = {m% + m% - m^J - {m\ + m\ - ml- J + 2ml (89) 



59 



can be as large as 0.6 x 10^ GeV^. This combination would vanish in the absence of the 
see-saw scale. There are also similar other mass combinations that have similar effects. 



2. Origin of matter. 

A very interesting aspect of the seesaw mechanism is the possibility that the heavy 
right handed neutrino decays and CP violation in the lepton sector may provide a way to 



understand the origin of matter - baryon asymmetry of the Universe 126 |. 

The original scenario consists of out-of-equilibrium CP violating decays of the RH neu- 
trinos N —>■ I + H which lead the to production of a leptonic asymmetry. This asymmetry 
is partly transformed to baryon asymmetry by sphaleron processes (which conserve B - L 
but violate B+L). 

One of the goals of this discussion is to learn about the right handed neutrinos and the 



nature of leptonic CP violation from the condition of successful leptogenesis |l27l | . Although 
there are many possible ways to achieve successful leptogenesis, e.g., resonant leptogenesis, 
non-thermal leptogenesis, etc., we restrict our discussion to the simplest case of a hierarchical 
pattern for right handed neutrino masses, i.e. Mi < M2,3 and thermal leptogenesis, and 
outline the consequences for neutrinos. 

The first implication of leptogenesis for the right handed neutrino spectrum comes im- 
mediately from the out-of-equilibrium condition for their decay: 

F, < H{M,) ^ v^Tj^, (90) 



where Fi is the decay rate 

^ {Y,Yj)uM, 



(91) 



»7r 

and Mi is the mass of the ith RH neutrino H{Mi) is the expansion rate of the Universe in 
the epoch with temperature T ~ Mj, and is the number of relativistic degrees of freedom 
in the epoch T. This condition leads to a lower bound on the mass of the RH neutrino 

M, > (92) 



One can get some idea about the required values of the masses Mi, e.g. assuming the up 
quark Yukawa couplings as a guideline for the Dirac neutrino couplings: Yj^^ik ~ Triujv^k- 

60 



Then the out-of-equihbrium conditions in eg. (1921) would imply that Mi > 10^ GeV, M2 > 
10^^ GeV and M3 > 10^^ GeV. So, the lepton asymmetry is assumed to be produced by the 
decay of the lightest A^i. 

In the leptogenesis scenario the baryon asymmetry, tib = n^/s, where Ub is the number 



density of baryons and s is the entropy density, can be written as [127 1 



Here rii is the number density of the RH neutrinos, ei is the lepton asymmetry produced in 
the decay of A^^i and ki is the wash out factor which describes the degree of out-of-equilibrium 
condition; the factor 8/23 is the fraction of the L— {B — L) asymmetry which is converted 
to the baryon asymmetry by sphalerons. The quantities rii, ei and /ti are all the functions 
of the RH neutrino masses and the Dirac type Yukawa couplings. So, the bou nds on th e 



127 . 



128 | 



neutrino parameters can be obtained from simultaneous analysis these quantities 
Indeed, some more information can be gained by analyzing the magnitude of the lepton 
asymmetry ei in terms of the Yukawa couplings Yy-. 

^1 = ^yyU E [iyX)lu] f , (94) 

where in the case of hierarchical mass spectrum of the RH neutrinos, x = M^/Ml ^ 1, 
and the function /(x) can be approximated as /(x) ~ —3/{2y/x) simplifying the above 
expression. Because rjB = 6.3 • 10^^°, and ki is roughly of order 10~^ (although it depends 
strongly on parameters in the model), we must have ei > 10~^, which puts according to 
a constraint on the flavor structure of Y^. 

Another constraint on Y^ follows from consideration of the decay rate of A^^i, which can 
be rewritten as fl9Tll 



r.^fi^?, . (95) 



This rate controls the initial abundance of A^^i and also the out-of-equilibrium condition. 
Successful leptogenesis restricts the values of rhi and Mi as shown in Fig. O According to 
results of fig. M this standard scenario implies a lower bound on the lightest RH neutrino 
mass Ml > 10*^ GeV and correspondingly gives an upper bound on the light neutrino 
masses thereby essentially excluding the degenerate spectrum for the case of a ty pe 



see-saw. These constraints can be avoided/weakened if one assumes type II seesaw [129 1 
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FIG. 9: Contour plot of the baryon to photon ratio produced in thermal leptogenesis, as a function 

128l |. The decay asymmetry ei was taken to be 10~^. The three (blue) 



of Ml and rhi, from 



close-together lines are the observed asymmetry. The horizontal contours, for small rhi assume a 
thermal A'^i abundance as the initial condition. 



and/or some specific flavor structures of the Yukawa couplings [l30|]. The bound can be also 
weakened in the case of strong degeneracy of RH neutrino masses Mi ^ M2 which leads to 



enhancement of the asymmetry ei (resonance leptogenesis 131 1]) and therefore allows for 
smaller ki. Consequently, the washout (out of equilibrium) conditions relaxes the bound on 
Ml. 



It was proposed recently that cosmological dens ity p erturbations can be generated by 
the inhomogeneous decay of right-handed neutrinos 132] • That requires coupling of the RH 
neutrinos with a scalar field whose fluctuations are created during inflation. 



3. Lepton flavor violation as tests of seesaw. 

Once one includes the right handed neutrinos in the standard model so that neutrinos 
acquire masses and mixings, lepton flavor changing effects such as /i — e + 7, r — ej, 
T ^ fi'-f appear. However, a simple estimate of the one-loop contribution to such effects 
leads to an unobservable branching ratio (of order ~ 10~^°). 
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The situation changes drastically as soon as the seesaw mechanism is embedded into 
the supersymmetric models. Flavor changing effects arise from the mixings among sleptons 
(superpartners of leptons) of different flavors caused by the renormalization group corre ction s 
which via loop diagrams lead to lepton flavor- violating (LFV) effects at low energies 133 1. 

The way this happens is as follows. In the simplest = 1 supergravity models, the 
supersymmetry breaking terms at the Planck scale are taken to have only few parameters: 
a universal scalar mass mo, universal A terms, and one gaugino mass mi/2 for all three 
types of gauginos. Clearly, a universal scalar mass implies that at Planck scale, there is 
no flavor violation anywhere except in the Yukawa couplings. However, as we extrapolate 
this theory to the weak scale, the flavor mixings in the Yukawa interactions induce flavor- 
violating scalar mass terms. In the absence of neutrino masses, the Yukawa matrices for 
leptons can be diagonalized so that there is no flavor violation in the lepton sector even 
after extrapolation down to the weak scale. However, when neutrino mixings are present, 
there is no basis where all leptonic flavor mixings can be made to disappear. In fact, in the 
most general case, of the three matrices: Yi - the charged lepton coupling matrix, Yi, - the 
RH neutrino Yukawa coupling and - the matrix characterizing the heavy RH neutrino 
mixing, only one can be diagonalized by an appropriate choice of basis and the flavor mixing 
in the other two remains. In a somewhat restricted case in which the right handed neutrinos 
do not have any interaction other than the Yukawa one and an interaction that generates 
the Majorana mass for the right handed neutrino, one can only diagonalize two out of the 
three matrices {i.e., Y^, Yi and Mr). Thus, there will always be lepton flavor violating terms 
in the basic Lagrangian, no matter what basis is chosen. These LFV terms can then induce 
mixings between the sleptons of different flavor and lead to LFV processes. The smallness 
of neutrino mass does not matter here. 

In the flavor basis, searches for LFV processes such as r — >■ /i -|- 7 and/or /i — >■ e + 7 can 
throw light on the RH neutrino mixings and/or family mixings in rriD, as has already been 
discussed. 

In the absence of CP violation there are at least six mixing angles (nine if m^i is not 
symmetric) in the seesaw formula, only three of which are observable in neutrino oscillation. 
To get useful information on the fundamental high scale theory from LFV processes, it is 
often assumed that is diagonal so that one has a direct correlation between the observed 
neutrino mixings and the fundamental high-scale parameters of the theory. The important 
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point is that the flavor mixings in then reflect themselves in the slepton mixings that 
lead to the LFV processes via the RGEs. 

To give a typical estimate of the magnitude of lepton flavor violation in seesaw models, 
we can make a simple ansatz of equal RH neutrino masses and assume CP conservation. 
The slepton mixing deflned by ^LL,ij = can be estimated from the renormalization 
group equation to be 

where M is the seesaw scale. Using M ~ 10^^ GeV or so, one flnds that the branching ratios 
for /i — »• e + 7 and r — >• /i + 7 depend on the slepton masses as and go down as slepton 
masses increase as can be seen from 

For the masses fh in the (200 - 500) GeV range, the e + 7 ) can be above 10~^^, a 



value that can be probed by the MEG experiment in progress 134j |. Similarly for the same 
slepton masses, R{t — > /i + 7) can be in the range of 10~^ — 10~^ or so. 

VII. NEUTRINO MASS AND NEW PHYSICS: "TOP-DOWN" 
A. Neutrino mass and Grand Unification. 

br p hysics beyond the Standard Model is supersymmetric grand 



One of the major ideas 

uniflcation (SUSY GUT) [135|. It is stimulated by a number of observations that are in 
accord with the general expectations from SUSY GUT's: (i) A solution to the gauge hier- 
archy problem, i.e. why ^ (h) uniflcation of electro-weak, i.e., SU(2)l x U{1)y 
and strong SU{?>)c gauge couplings assuming supersymmetry breaking masses are in the 
TeV range, as would be required by the solution to the gauge hierarchy problem; and (iii) a 
natural way to understand the origin of electroweak symmetry breaking. 

As noted above, the closeness of the gauge coupling uniflcation scale of approximately 10^^ 
GeV to the estimate of the seesaw scale from atmospheric neutrino data of M3 ~ 10^^ GeV 
suggests that the see-saw scale could be the GUT scale itself. Thus the smallness of neutrino 
mass goes quite well with the idea of supersymmetric grand uniflcation. However, in contrast 
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with items (i) through (iii) hsted above, the abundance of information for neutrinos makes 
seeing whether the neutrino mixings indeed fit well into simple SUSY GUTs a nontrivial 
exercise. 

The simplest GUT group is SU(5). Because the basic matter representations of SU(5), 
5 © 10, do not contain the right handed neutrino, one must extend the model by adding 
three right-handed neutrinos, one per generation. The problem then is that the Majorana 
mass of the gauge singlet right handed neutrino is unconstrained and can be the same as 
the Planck mass which will make it difficult to accommodate the neutrino data. The right 
handed neutrino mass fine-tuning question, i.e., of why Mr <^ Mpi arises again. However, 
if one includes the 15-dimensional Higgs boson, then the Yukawa interaction Em^m^^H in 
the superpotential leads to the coupling LLA^ where Aj, is the SU{2)l triplet in the 15- 
Higgs. The A° field acquires a VEV of order v'^^./XMu |136| . where A is a typical coupling 
parameter for the Higgs fields among themselves and Mjj is the scale of grand unification. 
This leads to neutrino masses (as in type II seesaw) of the right order of magnitude to 
explain the data. 

rn 

However, if one considers the SO (10) group [76], then its basic spinor representation 
automatically contains the right-handed neutrino along with the other fifteen fermions of 
the Standard Model (for each family). In order to give a mass to the right handed neutrino, 
one must therefore break SO(IO) symmetry (or more precisely, the B-L subgroup of SO(IO)). 
This naturally solves the right handed neutrino mass fine-tuning problem. Thus, one could 
argue that small neutrino masses have already chosen 5*0 (10) GUT as the most natural 
way to proceed beyond the Standard Model. Therefore 5*0(10) has rightly been the focus 
of many attempts to understand neutrino mixings. 

The 5*0(10) SUSY GUT models can be broadly classified into two classes. One class 
of models that employ the 16-Higgs representation to give mass to the right handed neu- 
trinos and another that employs a 126 Higgs. We outline below their major features and 
differences. 

As noted, one of the features that distinguishes S*O(10) from SU(5) is the presence of 
local B — L symmetry as a subgroup, and the two classes of the S*O(10) models mentioned 
above differ in the way the B — L symmetry is broken: Breaking by 16// Higgs field gives 
A{B — L) = 1, whereas 126// leads to A{B — L) = 2. In the first case the right-handed 
neutrino mass necessarily arises out of a nonrenormalizable coupling, whereas in the second 
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case it arises from a renormalizable coupling. In addition, the breaking of S — L by 16 Higgs 
leads necessarily to low energy MSSM with R-parity breaking so that the model cannot have 
cold dark matter without additional assumptions such as matter parity which forbids specific 
couplings such as (16m,)^16//, where 16^ stands for the matter spinor. 

However, 126// breaking of B — L preserves R-parity at low energies, so that the low- 
energy MSSM that derives from such an SO{10) has a natural dark matter candidate, i.e. 
the lightest SUSY particle. 

Because 5*0(10) contains the left-right symmetric group as a subgroup, it can have either 
a type II or type I seesaw formula for neutrino masses depending on the details of symmetry 
breaking and parameter ranges of the theory. For instance, in the 16h based models, the 
type II seesaw term is negligible and therefore the neutrino masses are dictated by type I 
seesaw formula. In contrast, in 126 Higgs models, the neutrino mass can be given by either 
the first term or the second term in the general seesaw formula, or both. 



1. A minimal 126-based SO(IO) model. 

Because 16 ® 16 = 10 © 120 © 126, the natural minimal model considers all the three 
Higgs fields and couples them to the matter 16. A simpler model contains only 10 © 126 
in which case there a re on ly two Yukawa coupling matrices: (i) h for the 10 Higgs and (ii) 



/ for the 126 Higgs |l37l |. SO(IO) has the property that the Yukawa couplings involving 
the 10 and 126 Higgs representations are symmetric. Therefore, to further reduce the 
number of free parameters, one may assume that Yukawa couplings are CP-conserving, and 
CP-violation arises from other sectors of the theory {e.g. squark masses). In a certain 
basis, one of these two sets of Yukawa coupling matrices is diagonal, and the Yukawa sector 
has only nine parameters. Noting that the (2,2,15) submultiplet of 126// as well as the 
(2,2,1) submultiplet of 10 h each have a pair of standard model doublets that contrib ute t o 



charged-fermion masses, one can write the quark and lepton mass matrices as follows 



137j : 



Mu = hn^ + fvu, Md = hKd + fvd, (98) 
Ml = hKd - 3fvd, rriD = hnu - ^fvu, (99) 



66 



where Ku^d are the VEV's of the up and down standard model type Higgs fields in the IOh 
multiplet and Vu,d are the corresponding VEV's for the same doublets in 126h- Note that 
there are 13 parameters (nine parameters in the Yukawa couplings noted above and four 
VEV's for the four MSSM doublets in the 10, and 126 Higgs fields) in the equations above, 
and 13 inputs (six quark masses three lepton masses and three quark mixing angles and the 
weak scale). Thus, all parameters of the model that go into fermion masses are determined. 

To generate the light neutrino masses, we use the seesaw formula in eq. flHOj) . where / is 
the same 126 h Yukawa coupling as ab ove. Thus all parameters that give neutrino mixings 



except an overall scale are determined 138 1. 



To see how large mixings arise in this model let us assume that the seesaw type II gives 
the dominant contribution to the neutrino mass, so that 

= fv,, (100) 

where is the VEV of EW triplet from 126//. Then, using eq. fllOOl) . we obtain 

= ^{Md-Mi). (101) 

Furthermore, note that the minimality of the Higgs content leads to the following sum rule 
among the mass matrices: 

kMi = rMd + Mu. (102) 

We then find using the known mass and mixing pattern for quarks, that 

fx' A3 
A3 A2 A2 



Md^i ^ nib- 



(103) 



VA3 A2 1 ; 

where A ~ sin^^c" = 0.22 and the matrix elements are supposed to give only the approximate 
order of magnitude. An important consequence of the relation between the charged lepton 
and the quark mass matrices in eq.f ll03p is that the charged lepton contribution to the 
neutrino mixing matrix, is f/; ~ 1 + 0(A), close to the identity matrix. As a result large 
neutrino mixings must arise predominantly from the neutrino mass matrix given by the type 
II seesaw formula. In the actual calculations of course the charged lepton mixings are also 
taken into account. The phenomenological fact that rrif, — 171^ ^ ttIt-X^ for a wide range of 
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values of tan/3 and Eqs. fllOip . (11030 imply that the neutrino mass matrix takes roughly the 
form 

/A3 A3 An 



=-^{Md-Mi)^ mo 



A3 A2 A2 

,3 ;)^2 \2 



(104) 



X' AV 

where, except for numbers of order one, the entire neutrino mass matrix is characterized 
by the Cabibbo angle alone. It is easy to see that both Ou (the solar angle) and 623 (the 



atmospheric angle) are now large 
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140|. 



The main point illustrated by this model is that the large neutrino mixings need not be 
a consequence of symmetries but rather could arise dynamically out of 6 — r unification at 
high scale. Note that this requires the choice of small /sa which was however made to fit 
the quark sector and not to "fix" the neutrino mixings. Of course, one must understand the 
flavor structure of the h and / Yukawa couplings {e.g. why /sa is so small) from a higher 
scale theory. 

There are various ways to incorporate CP violation into these models (see 
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142|). 



One could simply assume that the Yukawa couplings are complex. However, in this case the 
simple connection between b — r u ni ficat ion and large neutrino mixing is lost. In the pure 



type II seesaw SO (10) models jl39[ | 



140l | the fermion mass fits give the CKM phase in the 



second quadrant whereas in the mixed (type I plus typ e II) seesaw generalization, one can 
get a good fit to CP violation with correct CKM phase 1411 ]. 



Coming to the rest of the Higgs sector, one can define the minimal SO(IO) model |144| . |145 | 
as the one with the complete Higgs set 10 ©126 ©126 ©210 for all symmetry breakings and 
fermion masses. In this model it is hard to get the type II seesaw to dominate, a pr oblem 



easily cured by adding a 54 Higgs without affecting the fermion mass discussion 



146j. The 



parameter range of th e min imal model that fits fermion masses also has trouble giving rise 



to coupling unification 



147l | . Thus, it seems that the non-fermionic part of the Higgs sector 



of the 126 type models has to be extended to keep fermion mass fits. 

An alt erna t ive w ay to incorporate CP violation is to consider the Higgs set 10, 126 
and 120 



142 



143| . with all couplings constrained by CP invariance. Despite five extra 



parameters over the minimal model, one can retain the understanding of large atmospheric 
mixing as a consequence of 6 — r unification and a prediction for ^^13 as well as the Dirac 
phase of neutrinos in the observable range. 
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2. 16 Higgs-based SO(IO) models. 



The other class of SO(IO) models for neutrinos that ha s bee n wid e. 
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149 . 



150 



y d iscussed in 



151^. An ad- 



the literature includes just lOj^, 16/^-, IQh and 45^^ 
vantage of these models is that they use low dimensional Higgs multiplets. However, 
because the only renormalizable term in these models is 16ml6ml0//, they can nei- 
ther explain the observed quark and lepton masses nor can they explain the neutrino 
masses. One has to therefore include higher-dimensional operators in the Yukawa cou- 
pling such as 16ml6ml6Hl6H, 16„16ml6i/16H, IQrnXQ n^^mX^H, l^rnXQu^^mX^H , 

16ml6ml0i/45i^, where 16m stands for various fermion generations. Of these, the first two 
give symmetric Yukawa couplings, the next two have no symmetry property, and the last one 
can be both symmetric as well as antisymmetric. Because each coupling is a 3 x 3 matrix, 
there are many more free parameters in such models than observables. A strategy employed 
is to impose additional discrete symmetries to reduce the number of parameters. This and 
the fact that one can have large R-parity violation are drawbacks for these models. 

However, these models have certain advantages: (i) it is possible to implement the 
doublet-triplet splitting in a simple way such that the low energy theory below the GUT 
scale is the MSSM and (ii) the threshold corrections to the gauge couplings are not excessive, 
so that no particular constraint on symmetry breaking is necessary for the gauge couplings 
to remain perturbative. Another distinction from the 126iif-based models is that the type 
II seesaw contribution to neutrino masses is small in this model. The MSSM Higgs doublet 
fields, i.e. Hu,d, are linear combinations of the doublets in 10^, 16^ and 16//. The right 
handed neutrino mass arises from the 16ml6ml6//16// couplings when u'^ component of 



16// acquires a VEV. Typically, one uses the lopsided mechanism 



1201 ] to generate large 



atmospheric neutrino mixing. Most 16-Higgs based SO(IO) models lead to a small val ue fo r 



^13, although it is possible to have variations of the model that lead to a bigger value 



151 |. 



B. Grand unification and flavor symmetry 

Although the hypothesis of grand unification goes naturally with the seesaw scale, the 
detailed flavor pattern, i.e., the hierarchical mass and mixing among quarks and large mix- 
ings for leptons, is perhaps suggestive of some kind of flavor symmetry connecting different 
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generations. A possible symmetry group for such models that unify quark and lepton flavor 
textures while implementing the see saw mechanism coul d be, e. g. ^_S0 (10) family, where 
G family cau be eit her S O (3), SU(3) |152l |. or SU(2) |l54| or f/(l) |153l ] group, or a discrete 



group such as 5*4 



155| . Z2 or y44, all of which have been attempted. The groups such as 



SU{3), S0{3) as well as S4 and ^4 have an advantage over the U{1) and Z2 groups since 
they have 3 dimensional representations into which the three families can fit unlike the other 
groups. 

The main feature of these models is that in the case of abelian discrete group one can 
reproduce the flavor structure selecting the Yukawa couplings, whereas in the case of non- 
abelian ones, the problem shifts to VEV alignment and particular form of the scalar poten- 
tial. This generally requires large number of Higgs fields with specific couplings. However, 
this appears to be a straightforward and promising direction for both quark-lepton and flavor 
unification and better models must be pursued. 

It is also worth noting that if simpler models such as the minimal SO(IO) model with 
126 H discussed above are experimentally favored, we must find a natural way based on 
some higher symmetry to generate the necessary form of the 126// Yukawa coupling / as in 
eq.([IOl. 

Research along this line are mostly at an exploratory stage but it is probably fair to 
conjecture that such a unified approach is more likely to succeed if the neutrino mass hier- 
archy is established to be normal rather than inverted since the unification group connects 
all fermion textures. Furthermore, a strong signal of an underlying symmetry would be 
a degenerate spectrum. Examples of such symmetries which in conjunction with type II 
seesaw lead to degenerate spectrum have been discussed in the literature 
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156 |. 



C. Grand unification and proton decay 

Since grand unified theories connect quarks and leptons, most such theories predict an 
unstable proton and therefore one could use proton decay as a signal of the specific nature 
of the grand unified theory. In supersymmetric theories since the dominant contribution to 
proton decay arises from dimension five operators which involve Yukawa couplings respon- 
sible for flavor structure of fermions, one may also hope to learn about the fermion textures 
from the proton decay modes. 
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In the context of SO 
in both 16h based 



'10) models, the predictions for p roto n decay have been studied 



157( 1 as well as 126h based models 158l |. Both cases have ty pical 



predictions for distinguishing models, e.g. the cannonical p vK^ in the case of Ref. 157| 
and n —>■ n^u in the 126h case at an observable level. 

It must however be stressed that a true test of grand unification would be the discovery 
of the gauge mediated proton decay mode p e+7r° which is completely model indepen- 



dent. The present lower limit on the partial lifetime for this mode is 5 x 10'^^ years 159 |. 
For supersymmetric theories however they are expected to be at the level of 10'^^ years or 
longer and are therefore beyond the reach of experiments with conventional technology. The 
situation is more hopeful for nonsupersymmetric theories. For example in two step non- 

(2)t. X SU{2)r X SU{A)^ SM, the lifetime is 



SUSY SO(IO) theories where 50(10) ^ SU 
predicted to be 1.44 x io32 i±o-7±i o±i.9 ygg^rs 160| (where the uncertainties include threshold 
effects and gauge coupling uncertainties etc.). The intermediate scale in these theories is 
around 10^^'^ GeV which is of the right order to explain the neutrino masses via the seesaw 
mechanism. Other non-SUSY exa mple s are minimal SU(5) models with a 15-plet to give 



neutrino masses by type II seesaw 16 1| . In these theories predicted proton life time can be 
close to the present limit but in any case has an upper limit of 1.4 x 10^^ yrs. 

Clearly, proton decay search will be a crucial part of the search for grand unification that 
explain neutrino masses. 

D. Neutrinos and extra dimensions 

One of the important predictions of string theories is the existence of more than three 
space dimensions. For a long time, these extra dimensions were believed to be small and 
therefore practically inconsequential as far as low energy physics is concerned. However, 
recent progress in the understanding of the nonperturbative aspects of string theories has 
opened up the possibility that some of these extra dimensions could be large without con- 
tradicting observations. In particular, models in which some of the extra dimensions have 



sizes as large sub-millimeter and the string sca^ 



e is in the few TeV range have 



The basic assumption of these 



attracted a great deal of phenomenological attention 
models, inspired by the D-branes in string theories, is that space-time has a brane-bulk 
structure, where the brane is the familiar (3-|-l)-dimensional space-time, with the stan- 
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dard model particles and forces residing in it, and the bulk consists of all space dimensions 
where gravity and other possible gauge singlet particles live. One could of course envision 
(3+d+l)-dimensional D-branes where d space dimensions have miniscule (< TeV~^) size. 
The main interest in these models has been due to the fact that the low string scale provides 
an opportunity to test the models using existing collider facilities. 

In general the extra dimensional theories can be divided into three broad classes: (i) very 
small size fiat extra dimensions (r ~ M[J^ or so); (ii) large fiat extra dimensions (i.e. r ~ 
millimeter and (iii) warped extra dimensions of Randall-Sundrum type. 

In models with MjJ^ sized extra dimensions, one can implement the seesaw mechanism 
to generate neutrino masses. These models fit in very well with the conventional grand 
unified theories. These models have become popular as a way to providing an alternative 
resolution of the doublet-triplet splitting problem of grand unified theories via orbifold com- 



pactification 



163| . As far as the flavor problem goes, if all the flavors are in the same brane. 



the presence of extra dimension does not distinguish between them and therefore does not 
throw any light on this issue. There are however mod els w here different fermion genera- 
tions are put in different locations in extra dimensions 164| which then leads to nontrivial 
flavor structure and a possible way to approach the flavor problem. Usually however extra 
assumptions such as symmetries are needed to get realistic models. 

Coming to models with large extra D models (case (ii)), a major challenge to them 
comes from the neutrino sector. There are several problems: (i) how to understand the 
small neutrino masses in a natural manner since the see-saw mechanism does not work here 
due to lack of a high scale; (ii) second problem is that if one considers only the standard 
model group in the brane, operators such as LHLH/M^ could be induced by string theory 
in the low energy effective Lagrangian. For TeV scale strings this would obviously lead to 
unacceptable neutrino masses. 



One mechanism suggested in Ref. |165[ | is to postulate the existence of gauge singlet 
neutrinos, ub, in the bulk that couple to the lepton doublets in the brane and additionally 
demand the theory to be invariant under the B — L symmetry so that the higher-dimensional 
operator LHLH/M^ is absent. In four dimensions the Yukawa couplings and consequently 
Dirac masses turn out to be suppressed by the ratio M*/Mp£, where Mp^ is the Planck mass. 
The latter is now an effective parameter related to the volume of the extra dimensions, 
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Vd = {'2nYRi...Rd, and the fundamental scale as 

= M^+'^Vd. (105) 

This suppression is sufficient to explain small neutrino masses and owes its origin to the 
large bulk volume in comparison with width of the brane The volume suppresses 

the effective Yukawa couplings of the Kaluza-Klein (KK) modes of the bulk neutrino to the 
brane fields. 

The appearance of the suppression can be shown using one extra dimension with coordi- 
nate y and radius R. The full action involving the i>B{x,y) can be written as 

/h — 
d'^xdy[iuB7i^d''uB + ii'BL{x,y)dyUBR{x,y) + -j=5{y)LHuBR{x,y) + /i.e.], (106) 



where /i = 0, 1, 2, 3 and H denotes the standard model Higgs doublet. By expanding the 
bulk field in the Fourier series we obtain 

Z/p(x, V) = —p^ V^n^ COS (107) 

where v'^^ is the fcth KK mode and the prefactor follows from normalization of the wave 
function. Then, according to eq. fll06p the effective 4-dimensional Dirac coupling of the 
neutrino z/j^^ equals k = h/ \/2ttRM^. Generalization for the case of d extra dimensions is 
straightforward: 2ttRM^, —>■ VdMf. Now using the relation between the four and 4 + d - 
dimensional Planck masses in eq. (11051) we get k = h^^, independent of the number and 
configuration of the extra dimensions. After standard model gauge symmetry breaking, this 



generates a Dirac mass for the neutrino |165l | given by 



m = . 108 

Mpe ^ ' 

For M, ~ 10 - 100 TeV, eq. (fT08D leads to m ~ (10"^ - lO"^)/^ eV. Because /i is a five- 
dimensional coupling, its value could perhaps be chosen to be ~ 10 in which case we get 
neutrino mass in the range of interest in the discussion of neutrino oscillations. Furthermore, 
the usual left-handed neutrino is mixed with all the KK modes of the bulk neutrino, with 
the same mixing mass ~ \plm. Because the fcth KK mode has a mass = the 
mixing angle is given by \/2mR/k. Note that for R ~ 0.1 mm, this mixing angle could be 
relevant for the subdominant effects in solar and supernova neutrinos. 
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The above discussion can be extended in a very straightforward manner to the case 
of three generations. The simplest thing to do is to add three bulk neutrinos ascribing 
the generation label to all fermion fields. Now k becomes a 3 x 3 matrix. One can first 
diagonalize this by rotating both the bulk and the active neutrinos. The mixing matrix then 
becomes the neutrino mixing matrix U discussed in the text. After this diagonalization one 
can perform the KK expansion, which leads to mixing of the active neutrinos and the bulk 
towers. There are now three mixing parameters, one for each mass eigenstate denoted by 

= \/2miR, and the mixing angle for each mass eigenstate to the A;th KK mode of the 
corresponding bulk neutrino is given hy ^/k. 

In four dimensions the KK modes of the RH neutrinos show up as sterile neutrinos. 
The main feature is that there is an infinite number of such neutrinos with increasing mass 
and decreasing mixing. This can lead to peculiar effects in neutrino oscillations. Untill now, 
however, no effects are found which leads to the upper limits on and hence on the radius of 
the extra dimension i?, given a val ue o f the neutrino mass (or the coupling h). According 
to detailed analysis performed in 166[ | one has > 0.02 eV for a hierarchical, > 0.22 eV 
for an inverted, and > 4.1 eV for a degenerate neutrino spectrum. In generical, for all three 
cases the most stringent bound comes from the solar neutrino data 

Coming to the third type of extra D models i.e. the Randall-Sundrum scenario, where 
one invokes a warped extra space dimension, understanding small neutrino masses is less 
straightforward and has not yet reached a level where its detailed phenomeno logic al impli- 
cations can be discussed although some interesting attempts have been made |167 |. 

Essentially, theories of extra dimensions provide us with qualitatively new mechanism 
of generation of a small Dirac neutrino mass. There are different scenarios; however, their 
common feature can be termed the overlap suppression: the overlap of wave functions of the 
left-handed, i^Liv), and right-handed, i^R{y), components in extra dimensions (coordinate 
y). The suppression occurs owing to different localizations of the I'liv) and I'R^y) in the 
extra space. The effective Yukawa coupling is proportional to the overlap. Thus, in the 
large flat extra dimensional scenario described above ul is localized in the brane which has 
volume 1/{M^:Y in extra space, whereas ur propagates in the whole extra space volume Vn. 
Thus, the overlap equals the ratio of the two: {l/Mf)/{Vd) which is precisely the factor 
we have discussed above. In the Randall-Sundrum scenario, and ur are localized into 
two different branes and the overlap of their wave functions is exponentially suppressed. In 
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addition, extra dimensions can be the origin of light sterile neutrinos. 



VIII. BEYOND THREE NEUTRINOS: STERILE NEUTRINOS AND NEW 
PHYSICS 

An important part of our understanding of physics beyond the standard model involves a 
knowledge of whether there are only three light neutrinos t'e.^t.r or there are others. Known 
low-energy particle physics as well as cosmology constrain the number and properties of any 
extra neutrino. The fact that the measurement of the invisible Z-width at LEP is accounted 



for by three known neutrinos to a very high degree of accuracy [168| implies that any extra 
light neutrino must not couple to the Z-boson and hence not the W-boson either. Extra 
neutrinos are therefore termed sterile neutrinos (z/^)- 

Sterile neutrinos can communicate with the usual active particles via Yukawa interactions. 
Non-zero VEVs of the corresponding scalar bosons generate the Dirac-type mass terms which 
lead to mixing of active and sterile neutrinos. In turn this mixing may have important 
theoretical and phenomenological consequences. 

A. Phenomenology of sterile neutrinos 

The possible existence of sterile neutrinos and their mixing have interesting consequences 
in particle physics, astrophysics and cosmology. Most of the studies however give bounds 
on masses and mixing of these neutrinos (see Ref. 169| for a recent review). In particular, 
if mixing of sterile neutrinos with active neutrinos is strong enough, they can come into 
equilibrium in the early universe and affect the big bang nucleosynthe sis. Present data on 



primordial ^He, and ''Li abundances impose a constraint Neff < 1.5 170| on the effective 
number of sterile neutrinos, which were in equilibrium in the epoch of nucleosynthesis. This 
in turn leads to the bound on the active-sterile mixing ^5 as a function of the mass ms- 

Strong bounds on parameters of sterile neutrinos also come from structure formation in 
the Universe, from solar and supernova neutrinos, and from studies of the electromagnetic 
radiation in the Universe (because sterile neutrinos have a radiative decay mode). 

One may ask whether there is any need to introduce light sterile neutrinos. There are 
several reasons that are very suggestive: 
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i) Interpretation of the excess of e^n events observed in the LSND experiment 17l| in 
lations imply the existence of one or more extra sterile neutrinos with 



terms of z/^ — z/p, osci 



mass 1 — 5 eV 



172 



173| . Such an interpretation has its own problems. Furth ermore it 



174. Another 



contradicts result of analysis of the large scale structure (LSS) in the Universe 
possibility is decay of a relati vely heavy sterile neutrino with mass ~ 0.01 — 0.1 MeV 175 |. 
The MiniBooNE experiment 176[ | is testing the LSND result. 

(ii) Spherically asymmetric emission of sterile neutrinos with ma ss in the keV range during 
supernova collapses may explain the phenomenon of pulsar kicks 177 1. 



(iii) Sterile neutrinos with mass m« 
of the dark matter in the universe 178 
not compatible with this proposal 180 1 



■3 keV were proposed to be the warm component 



179(1 . However recent analysis of the LSS data is 



iv) Oscillations of sterile neutrino s in the early Universe can be the origin of the lepton 



asymmetry in the Universe 



179 



181 |. 



(v) Weak (statistically insignificant) indications of the presence of sterile states come 
from solar neutrino data: the low Homestake rate and absence of the upturn of the energy 
spectrum at low energies |l82[ |. 

If for some reason the existence of sterile neutrinos is confirmed, it will be a major 
revolution in the landscape of neutrino physics. We discuss some physics implications of 
this below, focusing mainly on the question of how to understand the lightness of sterile 
neutrinos in the context of extensions of the standard model. 



B. Sterile neutrinos and properties of active neutrinos 

Sterile neutrinos may have very small mixings for a given mass and therefore their astro- 
physical and cosmological effects may be unobservable. Despite of this, they can strongly 
influence the mass matrix of active neutrinos and therefore affect the implications of the 
established experimental results for fundamental theory jl83 |. 

Suppose the active neutrinos acquire {e.g., via the see-saw) the Majorana mass matrix 
TTia. Consider one sterile neutrino, 5*, with Majorana mass ms and mixing masses with 
active neutrinos mj^ = {mes,rn^s,^Ts) ■ If ^ '^as, then after decouphng of S the 
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mass matrix of active neutrinos becomes 

My = rua - 



ms 



(109) 



where the last term is the matrix induced by S. Let us consider some possible effects. 

The active-sterile mixing (induced matrix) can be the origin of large lepton mixing. 
Indeed, nta may have the usual hierarchical structure with small mixing. The mixing pa- 
rameters nias can be chosen in such a w ay that the resulting matrix M.^ leads to large or 
maximal mixing of active neutrinos 184l |. 

The induced matrix can be the origin of particular neutrino symmetries. Consider uni- 
versal couplings of a singlet field 5* with active neutrinos: m^g = mo(l, 1, 1) = 1712/ ^/S. Then 
the induced matrix has the form 



6m c 



m2 



(110) 



where D is the democratic matrix with all elements equal to one. Suppose that the original 
active neutrino mass matrix has the structure 

^0 ^ 
1-1. (Ill) 
-1 1 



V2 



\ 



j 



Then the sum = nia+Sms reproduces the mass matrix for the tri-bimaximal mixing 183 |. 
Second sterile neutrino can generate matrix eg. (11 111) . Clearly this changes the implications 
of the neutrino results, which would require existence of sterile neutrinos, flavor blindness 
of their couplings, etc.. Because S is outside the SM structure (with RH neutrinos) it may 
be easier to realize some particular symmetries for the induced matrix. 



C. On the origin of sterile neutrinos. 

Understanding the origin of light sterile neutrinos is a challenge. Their masses are not 
protected by the EW symmetry and some new physics (symmetries, dynamics) should exist 
to explain why they are light. 
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1. Mirror model for the sterile neutrino. 



An interesting scenario for physics beyond the standard model in which there is an iden- 
tical copy of both the forces and matter present side by side with known forces and matter 
has been discussed in literature. This new copy is called the mirror sector of the famil- 
iar universe. The mirror sector communicates with the familiar one only via gravitational 
interactions. This idea was originally proposed by 

Lee and Yang originally proposed this idea to maintain an exact parity symmetry in the 
full universe containii ig th e mirror sector even though in each sector parity is violated in 



its weak interactions 



185| . Such scenarios have emerged recently in the context of string 



theories, where one has -Eg x symmetry of matter and forces, with each acting on one 
10- dimensional brane world, and where under mirror parity one brane goes into another. 
They are completely consistent with what is known about the low energy particles and forces 
as well as the standard big bang model of the universe if one assumes that in the process 
of evolution of the universe, the reheat temperature of the mirror sector is somewhat lower 
than that of the visible sector. Many interesting phenomenological consequences can follow 
in generic versions of such a theory at low energies such as neutrino oscillations, the dark 
matter of the universe etc.. 

The mirror model was applied to the description of neutrino oscillation physics in 186| . 



where it was noted that if sterile neutrinos indicated by the LSND results are confirmed, 
one of the ways to explain their lightness is to postulate the existence of the mirror sector 
of the universe, in which case the mirror neutrinos can play the role of the sterile neutrinos 
and their lightness will follow from arguments similar to the familiar neutrinos, e.g., via 
mirror seesaw. Electroweak symmetry can generate their mixings via operators of the form 
LHL' H' / M ., where M could be the Planck mass realizing the possibility that the two sectors 
mix via gravitational interactions. In general, M could represent the mass of any standard 
model singlet particle. This model also has the potential to lead to sterile neutrinos in the 
keV range that mix with known neutrinos. 
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2. Other possibilities. 



Other models for the sterile neutrino include the possibility that itmay be a modulino 



- one of the standard model single t fie lds present in string models 
singlet fermions in the Eq models 



187|, one of the extra 



188| . or one of t he seesaw right handed neutrinos that 



becomes massless owing to leptonic symmetries 189l | such as Lg — L^ — Lr or fi — r exchange 
symmetry. A general feature of these models is that in the symmetry limit one of the SM 
singlet fermions remains massless and can be identified with the sterile neutrino. Its small 
mass and mixing with active neutrinos are generated via the terms that break the symmetry. 



IX. CONCLUSION 



Recent discoveries in neutrino physics have opened up a new vista of physics beyond 
the standard model. In this review we have attempted to provide a glimpse of what we 
have learned from the discoveries and what future experiments hold in terms how far this 
understanding can go. A broad theme is the appearance of new lepton flavor physics that 
was absent in the standard model with massless neutrinos and possibly important ramifica- 
tions for the flavor physics of quarks. The main areas we focused on were (i) understanding 
small neutrino masses; (ii) understanding the flavor structure of leptons that leads to large 
mixings and possible new symmetries implied by it, and (iii) some possible implications of 
the existence of new types of neutrinos. Of the several scenarios for understanding the small 
neutrino masses, the seesaw mechanism seems to have an advantage over others in many 
respects: (i) it provides a bridge to quark physics via grand unified theories; (ii) it gives a 
simple mechanism for understanding the origin of matter in the universe; and (iii) it has 
interesting low energy tests in the arena of lepton flavor violation and electric dipole moment 
of leptons. As far as the lepton flavor puzzle is concerned, although there are many interest- 
ing proposals, the flnal answer is far from clear and the next generation of experiments will 
very likely shed light on this issue. This process will quite possibly reveal new symmetries 
for leptons, which, in the broad framework of quark-lepton unification, may throw new light 
on the quark flavor structure. We summarized different ways to understand lepton mixings 
with and without the use of symmetries and discussed possible tests. Evidence for new 
neutrino species mixing with known neutrinos will be a new surprise in addition to the large 
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mixing surprise and will be another revolution. It could raise questions such as: Are there 
new quark species corresponding to the new neutrinos? What role do extra neutrinos play 
in the evolution of the universe, e.g., is there a mirror sector to the universe or are there 
extra dimensions? 

We emphasize that the field of neutrino physics is at an important crossroad in its 
evolution at the moment, and further advances will depend on how we answer the questions 
raised in this review. Some of the answers will very likely come from the proposed 
experiments that will test issues such as, is the neutrino its own antiparticle, how are the 
neutrino masses ordered, and what is the absolute scale of neutrino mass? Further precision 
measurements of neutrino parameters, as well as searches for new (sterile) neutrinos, are of 
fundamental importance. 

The work of R. N. M. is supported by the National Science Foundation grant number 
PHY-0354401. 
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